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This set of lectures is aimed at an audience who may be hearing the term 'anyon' 
for the first time. We shall start by explaining what the term 'anyon' means and why 
they are interesting. Just as fermions are spin 1/2, 3/2,...., particles obeying Fermi-Dirac 
statistics and bosons are spin 0,1,...., particles obeying Bose-Einstein statistics, 'anyons' are 
particles with 'any' spin obeying 'any' statistics. One may wonder why such particles have 
not been seen until now. After all, had they really existed in nature, they would have 
been just as familiar as the usual fermions and bosons. The explanation is quite simple. 
As will be seen in the course of the lectures, even theoretically, anyons can only exist in 
two space dimensions, whereas the real world is three dimensional. This naturally leads 
to the next question, 'why bother to study them at all ?'. The answer is that there do 
exist phenomena in our three dimensional world that are planar - systems where motion 
in the third dimension is essentially frozen. Anyons are relevant in the explanation of such 
phenomena. Besides, the study of anyons has led to a considerable improvement in our 
theoretical understanding of concepts like quantum statistics. 

The theoretical possibility of anyons was put forward as early as 1977 [1]. However, 
anyons shot into prominence and became a major field of research only in the last few years. 
One reason for this upsurge of interest was the discovery that the experimentally observed 
FQHE [2] had a natural explanation [3] [4] in terms of anyons. An even more dramatic rise 
in its popularity occured when it was discovered that a gas of anyons superconducts [5] [6], 
when it is coupled to electromagnetism. In fact, for a while, 'anyon superconductivity' was 
one of the top "candidate" theories to explain high T c superconductors [7]. Now, due to 
lack of experimental confirmation, the theory is no longer a 'hot' candidate, but interest in 
the field of anyons remains as high as ever. The rest of these lectures will involve a more 
detailed elaboration on the theme of anyons. 

We shall first start in Sec.(l) with basic notions of spin and statistics and understand 
why anyons can only exist in two spatial dimensions [1] [8]. Then we shall study a simple 
physical model of an anyon that incorporates fractional spin and statistics [9]. In Sec. (2), 
using this model, we shall solve some simple two anyon quantum mechanics problems, and 
see that the anyon energy eigenvalues actually interpolate between fermionic and bosonic 
eigenvalues. We shall also discover that even non-interacting two anyon states are not 
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simple products of single anyon states [1]. This is the crux of the problem in handling many 
anyon systems. Here, we shall concentrate on two approximations in which the many anyon 
problem has been tackled. The quantum statistical mechanics of a many anyon system has 
been studied [10] via the virial expansion of the equation of state, which is valid in the 
high temperature, low density limit. In Sec. (3), we shall first briefly review the classical 
and quantum cluster expansions and the derivation of the virial coefficients in terms of the 
cluster integrals. Then, using the results of the two anyon problems discussed in Sec. (2), 
we shall derive the second virial coefficient of the anyon gas [11]. The other approximation 
in which the many anyon system has been studied is the mean field approach which is valid 
in the high density, low temperature limit. In this approximation, every anyon sees an 
'average' or 'mean' field due to the presence of all the other anyons. Thus, the many body 
problem is reduced to the problem of a single particle moving in an 'average' potential. It 
is in this mean field approach that anyon superconductivity has been established [5] [6]. In 
Sec. (4), we shall study the mean field approach and derive anyon superconductivity in a 
heuristic way. 

Finally, we shall briefly touch upon two slightly more advanced topics, just to give 
a flavour as to why the study of anyons form such an interesting and relevant field of 
research today. The more formal topic deals with the formulation of a field theory of 
anyons. In Sec. (5), we shall introduce anyons in a field theory formalism using a Chern- 
Simons construction [10] and study an explicit example [12] of a Lagrangian field theory 
whose topological excitations are anyons. The second topic deals with the more physical 
question of applicability of anyon physics to condensed matter systems. In Sec. (6), we 
shall briefly indicate how anyons arise in the FQHE, which is a system of two dimensional 
electrons at low temperatures and in strong magnetic fields and show how the idea of 
statistics transmutation is used in novel explanations [13] [14] of the effect. We wish to 
emphasise here that realistic anyons do not exist in vacuum as has been assumed in the 
earlier sections, but actually arise as quasiparticles in a real medium. 
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1. Basic Concepts in Anyon Physics 



Let us start with spin in the familiar three dimensional world. We know that spin is 
an intrinsic angular momentum quantum number that labels different particles. The three 
spatial components of the spin obey the commutation relations given by 

[Si, Sj] = ieijkSk- (1.1) 

We shall show that these commutation relations constrain S to be either integer or half- 
integer. Let \s, m) be the state with S 2 \s, m) = s(s + l)|s, to) and S%\s, m) = m\s, m). By 
applying the raising operator S + , we may create the state 

S+\s,m) = [s(s + 1) -m(m + l)] 1/2 \s,m + 1) = \s,m'). (1.2) 

Requiring this state to have positive norm leads to 

s(s + l) -to(to + 1) > V to, (1.3) 

which in turn leads to 

m < s V m. (1.4) 
Thus, it is clear that for some value of m' = m + integer, m! > s unless s — m! - i.e., 

s — m = integer. (1.5) 

Similarly, by insisting that S~ = \s,m) have positive norm, we get 

s(s + 1) — m{m — 1) > V m, (1.6) 

which in turn implies that 

m> —s V to. (1.7) 

Once again, we construct the states, S~\s,m), (S~) 2 \s,m),... and to avoid to < — s, we 



have to set 



m — (— s) = integer. 



(1.8) 



Adding equations (1.5) and (1.8), we get 



2s = integer 



integer 



(1.9) 



s = 



2 



Thus, just from the commutation relations, we have proven that particles in 3+1 dimensions 
have either integral or half-integral spin. 

In two spatial dimensions, however, there exists only one axis of rotation (the axis 
perpendicular to the plane). Hence, here spin only refers to S3, which has no commutation 
relations to satisfy. For a given magnitude of S3, it can only be either positive or negative 
depending upon the handedness of the rotation in the plane. Since there are no commuta- 
tion relations to satisfy, there is no constraint on S3 and hence, we can have 'any' spin in 
two dimensions. For completeness, note that in one spatial dimension, there is no axis for 
rotation and hence, no notion of spin. 

The term statistics refers to the phase picked up by a wavefunction when two identical 
particles are exchanged. However, this definition is slightly ambiguous. Does statistics 
refer to the phase picked up by the wavefunction when all the quantum numbers of the 
two particles are exchanged (i.e., under permutation of the particles ) or the actual phase 
that arises when two particles are adiabatically transported giving rise to the exchange? In 
three dimension, both these definitions are equivalent, but not so in two dimensions. We 
shall concentrate on the second definition which is more crucial to physics and return to 
the first definition later. 

Let us first consider the statistics of two identical particles moving in three space 
dimensions [1] [8]. The configuration space is given by the set of pairs of position vectors 
(ri,r2). The indistinguishability of identical particles implies the identification (ri,r2) ~ 
(r2,i"i) - i.e., we cannot say whether the first particle is at ri and the second particle at 
T2 or the other way around. We shall also impose the hard-core constraint, ri 7^ T2, to 
prevent intersecting trajectories so that we can determine whether or not two particles have 
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been exchanged. However, as we shall see later, this constraint is unnecessary because for 
all particles (except bosons), there is an automatic angular momentum barrier, preventing 
two particles from intersecting, whereas for bosons, we need not know whether the particles 
have been exchanged or not, since the phase is one anyway. For convenience in constructing 
the configuration space, let us define the centre of mass (CM) and relative coordinates - 
R = (ri + T2)/2 and r = ri — Y2- In terms of these coordinates, the configuration space is 
(R, r) with r^O, and with r being identified with — r. This can be written as 



ftSl^). (1.10) 



Here R3 denotes the three dimensional Euclidean space spanned by R. The notation 
(i?3 — origin) for r implies that the origin r = is being dropped. Z2 is just the multiplicative 
group of the two numbers 1 and -1. Hence, dividing (i?3 — origin) by Z2 implies the 
identification of every position vector r in the relative space with its negative — r. To study 
the phase picked up by the wavefunction of a particle as it moves around the other particle, 
we need to classify all possible closed paths in the configuration space. Notice that the 
CM motion just shifts the positions of the two particles together and is independent of any 
possible phase under exchange. Hence, we are only required to classify closed paths in 



(gg ~ origin) = s 
Z2 



Instead of dealing with paths in S, for ease of visualisation, let us construct a simpler 
configuration space by keeping the magnitude of r fixed , so that the tip of r defines the 
surface of a sphere. Furthermore, the identification of r with — r implies that diametrically 
opposite points on the sphere are identified. Thus, configuration space is the surface of a 
sphere with opposite points identified as shown in Figs.(la,lb,lc). 
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(a) (b) (C ) 

Fig. 1 

Now, since we have eliminated coincident points, the wavefunction is non-singular and 
well-defined at all points in the configuration space. In particular, it is non-singular on the 
surface of the sphere. Hence, the phase picked up by the wavefunction under an adiabatic 
exchange of the two particles is also well-defined and does not change under continuous 
deformations of the path. Let us consider the possible phases of the wavefunction when 
the motion of the particles is along each of the three paths - A (no exchange), B (single 
exchange) and C (two exchanges) - depicted in Figs. (la, lb, lc). Path A defines a motion 
of the particles which does not involve any exchange. It is clearly a closed path on the 
surface of the sphere and can be continuously shrunk to a point. So this path cannot 
impart any phase to the wavefunction. Path B, on the other hand, involves the exchange 
of two particles and goes from a point on the sphere to its diametrically opposite point 
- again a closed path. Since the two endpoints are fixed, by no continuous process can 
this path be shrunk to a point. Hence, this path can cause a non-trivial phase in the 
wavefunction. However, path C which involves two exchanges, forms a closed path on 
the surface of the sphere, which, by imagining the path to be a (physical!) string looped 
around an orange (surface of a sphere), can be continuously shrunk to a point. So once 
again, the wavefunction cannot pick up any phase under two adiabatic exchanges. This 
leads us to conclude that there are only two classes of closed paths that are possible in 
this configuration space - single exchange or no exchange. Let rj be the phase picked up by 
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any particle under single exchange. The fact that two exchanges are equal to no exchanges 
implies that 



rf = +1 =>■ 7] = ±1. 



(1.12) 



Hence, the only statistics that are possible in three space dimensions are Bose statistics 
and Fermi statistics. 

Why does this argument break down in two space dimensions? Here, configuration 
space is given by 

(R 2 - origin) 



R 2 



Z 2 



1.13) 



where i?2 is two dimensional Euclidean space. Just as before, we ignore the CM motion 
and fix the magnitude of the relative separation, so that the configuration space can be 
visualised as a circle with diametrically opposite points identified (see Figs. (2a, 2b, 2c)). 
A 




(a) 



(b) 



(c) 



Fig. 2 

Here, however, several closed paths are possible. The path A that involves no exchanges 
can obviously shrunk to a point since it only involves motion along any segment of the 
circle and back. However, the path B that exchanges the two particles is just as obviously 
not contractible since the end-points are fixed. But even the path C, where both the 
dashed and the solid lines are followed in the clockwise direction (or equivalently both in 
the anti-clockwise direction) cannot be contracted to a point. This is easily understood 
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by visualising the paths as physical strings looping around a cylinder (a circle in a plane). 
Thus, if rj is the phase under single exchange, r/ 2 is the phase under two exchanges, rf is 
the phase under three exchanges and so on. All we can say is that since the modulus of the 
wavefunction remains unchanged under exchange, 77 is a phase and can be written as e t0 , 
where 9 is called the statistics parameter. This explains why we can have 'any' statistics 
in two spatial dimensions. 

The crux of the distinction between configuration spaces in two and three dimensions is 
that the removal of the origin in two dimensional space makes the space multiply connected, 
whereas three dimensional space remains singly connected. Hence, in two dimensions, it is 
possible to define paths that wind around the origin. This is not possible in three dimen- 
sions. Mathematically, this distinction is expressed in terms of the first homotopy group 
111, which is the group formed by inequivalent paths (paths that are not deformable to one 
another), passing through a given point in configuration space, with group multiplication 
being defined as traversing paths in succession and group inverse as traversing a path in 
the opposite direction. Thus, in two dimensions, 

ITi(2dim.config.space) = U t ^ 2 ~ ongm ) = U^RP^ = Z (1.14) 

Z2 

where Z is the group of integers under addition and RP\ stands for real projective one di- 
mensional space and is just the notation for the circumference of a circle with diametrically 
opposite points identified. The equality in the above equation stands for isomorphism of 
the groups so that the homotopy group of configuration space is isomorphic to the group 
of integers under addition. The different paths are labelled by integer winding numbers, 
so that the phases developed by the wavefunction are of the form, r] n , with n an integer , 
which in turn leads to the possibility of 'any' statistics in two dimensions. In contrast, in 
three dimensions, we have 

rii(3dim.config.space) = III ^ 3 ~ origm ) = Ui(RP 2 ) = Z 2 . (1.15) 

Z2 

Here RP2 stands for real two dimensional projective space and is the notation for the 
surface of a sphere with diametrically opposite points identified. Since, Zi has only two 
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elements, there exist only two classes of independent paths and thus, only two possible 
phases - fermionic or bosonic - in three spatial dimensions. 

The distinction between the phase of the wavefunction under exchange of quantum 
numbers and the phase obtained after adiabatic transport of particles is also now clear. 
Under the former definition, the phase rf after two exchanges is always unity, so that 
i] = ±1, whereas the phase under the latter definition has many more possibilities in two 
dimensions. Mathematically, the distinction is that the first definition classifies particles 
under the permutation group P/v, whereas the second definition classifies particles under 
the braid group Bjy. The permutation group (-P/v) is the group formed by all possible 
permutations of N objects with group multiplication defined as successive permutations 
and group inverse defined as undoing the permutation. It is clear that the square of any 
permutation is just unity, since permuting two objects twice brings the system back to the 
original configuration. Thus, particles that transform as representations of P/v can only be 
fermions or bosons. The braid group Bjy, on the other hand, is the group of inequivalent 
paths (or trajectories) that occur when adiabatically transporting N objects. For example, 
the trajectories shown in Fig. (3) 




(a) (b) ( C ) 

Fig. 3 

are elements of B4, because all of them are possible paths involving four particles. Fig. (3b) 
represents the identity element where none of the trajectories cross each other. Group 
multiplication is defined as following one trajectory by another as depicted in Fig. (4) 
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© 



Fig. 4 



and group inverse is defined as a reverse crossing, (Fig. (5)) 
1 — 



Fig. 5 



so that the product of a trajectory and its inverse leads to the identity as shown in Fig. (6). 





Fig. 6 



11 



Here, it is pictorially clear (see Fig. (7)), 




Fig. 7 

that the square, or indeed, any power of the trajectory representing the adiabatic exchange 
of two particles is not 1. Hence, particles that transform as representations of the braid 
group are allowed to pick up 'any' phase under adiabatic exchange. For completeness, we 
mention that more abstractly, the braid group Bn is defined as the group whose elements 
(trajectories) satisfy the following two relations depicted pictorially in Fig. (8) 




Fig. 8 

and Fig. (9). 
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Fig. 9 



The second relation is called the Yang-Baxter relation. It is clear that the braid group is 
a much richer group than the permutation group and leads to a much finer classification. 
For example, consider the trajectories a) and b) in Fig. (10). 

12 12 




1 2 

(a) 



Fig. 10 



1 2 
(b) 



The permutation group cannot distinguish between these two trajectories. In both cases, 
there is no permutation and hence the trajectory belongs to the identity representation, 
whereas the two trajectories are distinct elements of the braid group. 

Finally notice that in one spatial dimension, two particles cannot be adiabatically 
exchanged without passing through each other -i.e., without interacting. Hence any theory 
can be written in terms of bosons or fermions with appropriate interactions and there is no 
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real concept of statistics in one dimension. 

After all this abstract discussion, let us construct a simple physical model of an 
anyon [9]. Imagine a spinless particle of charge q orbiting around a thin solenoid along 
the z-axis, at a distance r as shown in Fig. (11). 




Fig. 11 



When there is no current flowing through the solenoid, the orbital angular momentum of 
the charged particle is quantised as an integer -i.e., 



l z = integer. 



;i.i6) 



When a current is turned on, the particle feels an electric field that can easily be computed 
using 

y (V x E)d 2 r ■■ 
where is the total flux through the solenoid. Hence, 



(1.17) 



E-dl = 2n\r\E e 



leading to 



E = 



2tt r 



[z x r). 



(1.18) 



(1.19) 



Thus, the angular momentum of the charged particle changes, with the rate of change being 
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proportional to the torque r x F -i.e., 



Hence, 



l z = r x F = r x qE = (1.20) 

Its 



= -I (1.21) 



is the change in angular momentum due to the flux <fi through the solenoid. In the limit 
where the solenoid becomes extremely narrow and the distance between the solenoid and 
the charged particle is shrunk to zero, the system may be considered as a single composite 
object - a charged particle-flux tube composite. Furthermore, for a planar system, there 
can be no extension in the z-direction. Hence, imagine shrinking the solenoid along the 
^-direction also to a point. The composite object is now pointlike, has fractional angular 
momentum and perhaps can be identified as a model for an anyon. However, this is too 
naive a picture. As we shall see later, in an anyon, the charge and the flux that it carries 
are related. So it is not quite right to think of the anyon as an independent charge orbiting 
around an independent flux. The charge is actually being switched on at the same time 
that the flux in the solenoid is being switched on. Hence, Eq.(1.20) needs to be modified 
to read 

4 = (1.22) 
Moreover, since q(t) = c(f>(t), for some constant c, we get 

Al z = - C 4-= & (1.23) 

so that the angular momentum of a charge-flux composite, with charge proportional to the 
flux, is less than what we originally computed. This is not surprising since the original 
computation overestimated the charge by keeping it fixed. Hence, our final physical model 
of an anyon is that of a charge-flux composite with charge q and flux <fi being proportional 
to each other and with a spin given by q<p/47r. In the next section, we shall see that this 
model also exhibits fractional statistics and complete its identification as an anyon. 
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Problems 



1. Fermions and bosons are one-dimensional representations of the permutation group 
Pjv and anyons are one-dimensional representations of the braid group Bjy. Can 
one have higher dimensional representations of P/v and P/v? (The answer is yes.) 
Think about the quantum mechanics and statistical mechanics of particles in these 
representations and see how far you get. 
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2. Quantum Mechanics of Two Anyon Systems 



In Sec.(l), we constructed a composite object which consisted of a spinless (bosonic) 
charge orbiting around a (bosonic) flux, and showed that when the charge is proportional 
to the flux, this object had fractional spin s = q(f)/4ir. To determine its statistics, we need 
to study the quantum mechanical system of two such objects. The Hamiltonian for the 
system is given by 

H _ (pi- gai) 2 ( P2 - ga 2 ) 2 

2m 2m 1 ' 1 



with 



(f) z x (ri - r 2 ) 

ai = 



2n I 1 " 1 " 1 " 2 ! (2 2) 

z x (r 2 - ri) 
and a 2 = , 

2lX \V\ — Y2\ 

where ai and a 2 are the vector potentials at the positions of the composites 1 and 2, due to 
the fluxes in the composites 2 and 1 respectively. Let us work in the centre of mass (CM) 
and relative (rel) coordinates - i.e., we define respectively, 

ri + r 2 Pi - P2 /(-j q\ 
R= — - — ^P = pi + p 2 , and r = n - r 2 => p = . (2.3) 

In terms of these coordinates, the Hamiltonian can be recast as 

H = + (P " gare ° 2 (2.4) 
Am m 

with 

z x r 

arel = 2^TF- (2 - 5) 

Thus, the CM motion, which translates both the particles rigidly and is independent of 
statistics, is free. The relative motion, on the other hand, which is sensitive to whether the 
composites are bosons, fermions or anyons, has reduced to the system of a single charged 
particle of mass m/2 orbiting around a flux <p at a distance r. Since the composite has 
been formed of a bosonic charge orbiting around a bosonic flux, the wavefunction of the 
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two composite system is symmetric under exchange and the boundary condition is given 
by 

VVelM + TT) =VrelM) (2.6) 

where ip Te \ is the wavefunction of the relative piece of the Hamiltonian in Eq.(2.4) and 
r = (r, 9) in cylindrical coordinates. 

Now, let us perform a (singular) gauge transformation so that 

a rel — > 4 el = a rel - VA(r, 9) (2.7) 

where A(r, 9) = ^9. This gauge transformation is singular because 9 is a periodic angular 
coordinate with period 9 and is not single valued. In the primed gauge, 

1 (f)9 

a reW = a reW ~ 'ggi^) = alld fl relr = a relr = (2.8) 

-i.e., the gauge potential completely vanishes. Hence, the Hamiltonian reduces to 

p2 2 

H = T~ + — (2-9) 
4m m 

which is just the Hamiltonian of two free particles. However, in the primed gauge, the 
wavefunction of the relative Hamiltonian has also changed. It is now given by 

V4i(r, 9) = e-^ rel (r, 9) = e~^% el (r- 9) (2.10) 

which is no longer symmetric under r — > — r since 

< el (r, 9 + tt) = e-^/Vki(r, fl) = e- iQ ^ el (r, 6»). (2.11) 

Thus, two interacting bosonic charge-flux composites are equivalent to two free particles 
whose wavefunctions develop a phase e - * 9 ^/ 2 under exchange - i.e., they obey fractional 
statistics. This completes the identification of charge-flux composites as anyons. Notice 
that the statistics phase a = q<f>/2 is in accordance with the spin of the composite, which we 
had earlier determined to be q<f)/4ir, so that the generalised spin-statistics theorem which 
relates the statistics factor a to the spin j = a/2n, is satisfied. 
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We have thus shown that anyons can either be considered as free particles with frac- 
tional spin obeying fractional statistics, or as interacting charge-flux composites, again with 
fractional spin, but obeying Bose statistics. The free particle representation is called the 
anyon gauge, and the interacting particle representation is called the boson gauge. 

Let us now solve the quantum mechanical problem of two free anyons. The Hamiltonian 
in Eq.(2.9) is simple, but the boundary condition in Eq.(2.11) is non-trivial. However, we 
have just seen that this problem is equivalent to the problem of two interacting charge-flux 
composites, described by the Hamiltonian in Eq. (2.4) with bosonic boundary conditions. 
The CM motion is trivial and directly solved to yield the energy eigenvalues and eigen 
functions as 

p2 

EcM = 4m" ^CM = e lP ' R . (2.12) 

To solve for the relative motion, we work in cylindrical coordinates, wherein the relative 
part of the Hamiltonian equation takes the form 

["-(It + ~h + - ^) 2 ]^el(r, 9) = E rel A el (r- 9) (2.13) 

m Or 1 r Or mr z 09 2n 

with the boundary conditions in Eq.(2.6). Since the Hamiltonian is separable in r and 9, 
the wave function factorises as ip Te \(r, 9) = 1Z(r)Yi(9). Hence, the angular equation reduces 
to 

^9't^ )2Ym = XYm (2 ' 14) 
whose solution for Yi(9), consistent with the boundary condition is given by 

Y l (9) = e ae , l = even integer, (2.15) 

which, in turn, leads to 

A = (/ + ^) 2 , 1 = even integer. (2.16) 
2% 

When the A eigenvalue is substituted in the radial equation, we get 

[--(1^ + -#) + A( Z + ?) 2 Mr) = E iel K(r). (2.17) 
m or 1 r or mr 1 2ir 

Notice that the net effect of the fluxtubes has been to add a factor q<p/27r to the angular 
momentum term, thus adding to the centrifugal barrier. This justifies our earlier state- 
ment in Sec.(l) that all particles, except bosons, have a centrifugal barrier which prevents 
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intersecting trajectories. Now, by defining mE re [ — k , we may rewrite the radial equation 
in the form 

which is easily identified as the Bessel equation with the solution 

H(r) = J ll+ ^(kr). (2.19) 

Hence the relative wavefunction of the composites, in the bosonic gauge, is given by 

V'relO", 0) = e M J^ l+ ^{kr), I = eveninteger. (2.20) 

The wavefunction for two anyons can also be expressed in the anyon gauge by performing 
a gauge transformation so that 

^ el = e^ l+ ^ e J ]l+ ^(kr), 1 = even integer. (2.21) 

This wavefunction is obviously anyonic, since it picks up a phase q<f>/2 under 6 — > 9 + n 
(i.e., under exchange of the two particles). 

Including the CM motion, the two anyon wavefunction can be written as 

^(R,r) = VcMWrelW = e ip - R e^ l +^ e J ll+ ^(kr) (2.22) 

or equivalently as 

^(R,r) = e iLQ J L {KR)e i{l+q ^ e J {l+q ^{kr). (2.23) 

In Eq.(2.23), we have expressed the CM motion also in terms of cylindrical quantum num- 
bers - L is the CM angular momentum and K = AmEcM labels the CM energy. The two 
particle wavefunction can be recast in terms of the original single particle coordinates ri 
and Y2 - i.e., ip(R, r) — ► ip(r\, r2). However, unless |^ is either integral or half-integral, the 
two particle wavefunction cannot be factorised into a product of two suitable one particle 
wavefunctions — 1*2) 7^ -? /'i( r i)'?/'2( r 2)- Since handling Bessel functions is inconvenient, 
this property will be demonstrated more explicitly in the next example of two anyons in a 
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harmonic oscillator potential. Furthermore, in principle, we should be able to prove from 
this example that the energy levels of a two anyon system cannot be obtained as sums of 
one anyon energy levels. This, again, is easier to see with discrete energy levels and will be 
explicitly proven in the next example. 

The second example that we shall solve explicitly is the problem of two anyons in a 
harmonic oscillator potential [1] with the Hamiltonian given by 



As in the case of two free anyons, the problem is separable in the CM and relative coordi- 
nates, in terms of which, the Hamiltonian can be rewritten as 



The CM motion is clearly independent of the statistics of the particles and involves just the 
usual quantum mechanical problem of a single particle of mass 2m in a two dimensional 
harmonic oscillator potential. The energy levels and wavefunctions are obviously well- 
known. However, we shall briefly recollect the familiar steps here, just to set the field for 
the relative motion problem which is sensitive to the statistics of the particles. For the CM 
motion, we work in the cylindrical (R, 6) coordinates and write 




P9 1 99 1 99 
+ — + -mo/rf + -muj z ri. 
2m 2 1 2 2 



(2.24) 




P , 2-d2 1 1 2 2 

h mu R + -muj r . 



m 4 



(2.25) 



H CM ^ CM (R,Q) = [ 



I , d 2 Id 




_L cP_ 

AmR 2 dQ 2 



+ mu 2 R 2 }^cM{R,Q)- (2.26) 



Am y dR 2 RdR 



The wavefunction factorises in R and and can be written as 



iJ>Cm(R,Q) = e lL »K CM (R) 



(2.27) 



with L being an integer for a single valued wavefunction. Defining K 2 
radial eigenvalue equation becomes 



= 4mEcM, the 



(2.28) 



21 



with a solution of the form 

oo 

n CM (R) = e- m " R2 J2anR n+s . (2.29) 

n—o 

By substituting this series solution in Eq.(2.28), we find that s = \L\ and that 

a n+ 2 K 2 - Amuj(n + \L\ + 1) 
= (n + 2) 2 + 2|L|(n + 2) ' 

so that n = even integer. Requiring the series to terminate leads to 

E CM = cu(n+ \L\ + 1) (2.31) 

which is the well known answer, since n + \L\ = p is an integer, and the modulus factor 
gives the appropriate degeneracies. 

The Hamiltonian for the relative motion, which does depend on the statistics of the 
particles, can also be solved in the same way. In the boson gauge, we have 

#rel^rel = [ ~ h 7™^ 2 ]^rel = ^rel^rel (2-32) 

m 4 

with a re i = (O,0/27rr) and ip re \(r, 9 + ir) = ip Te i(r,9). The wavefunction, once again, fac- 
tor ises in r and 6 and is given by 

Ael(r,9) = e ue K iel (r) (2.33) 

However, due to the boundary condition on ip Te \, I now has to be an even integer, analo- 
gous to the relative angular momentum quantum number for two free anyons . With the 
definition k 2 = mE, the radial equation is 

l^ + l§- r - ±(l + a/,) 2 - \mW + k 2 \KUr) = 0, (2.34) 
where we have substituted a = qcp/2. As in the case of two free anyons, we see that the net 
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(2.30) 



effect of statistics is to add a term to the centrifugal barrier. A series solution of the form 



oo 

„2 



ft rel ( r ) = e~ mu}r / 4 b n r n+s (2.35) 



n=0 



can be found to Eq.(2.34), with s — \l + a/ir\ and 

b n +2 k 2 — muj(n + \l + a/ix\ + 1) 
= (n + 2)2 + 2|/ + a/vr|(n + 2)' 



(2.36) 



so that n has to be an even integer. As before, the requirement that the series has to 
terminate leads to the energy levels given by 

E ve i = u(n+\l + a/n\ + l). (2.37) 

Notice that a = and a = ir give the usual energy levels for bosons and fermions respec- 
tively. The first few energy levels and their degeneracies are explicitly given by 



Eq = 


(1 + a/n)uj, 


deg = 


1, 


Ei = 


(3 — a/n)uj, 


deg = 


1, 


Ei = 


(3 + a/7i)u, 


deg = 


2, 


Es = 


(5 — a/n)uj, 


deg = 


2, 


E 4 = 


(5 + a/n)u, 


deg = 


3. 



(2.38) 



Thus, it is clear that the energy levels can be written as 

Ej = (2j + 1 + a/n)u, deg =j + l 
and Ej = (2j + 1 — a/n)u, deg = j, 



(2.39) 



where j is an integer. The energy levels can be plotted as a function of a as shown in 
Fig.(12). 
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06 = 



1 



2 



a = 7i 



Fig. 12 



Here a = corresponds to bosons and a = ir corresponds to fermions. The degeneracies 
of the levels are mentioned within brackets. From the figure, it is clear that the energies 
of the 2 particle system are a monotonic function of the a and change continuously as a 
changes from bosonic to fermionic value. Moreover, except at a — 0,7r/2,7r, the energy 
levels are not equally spaced and even at a = n/2, the spacing between energy levels is 
half that for bosons and fermions. We also note that the energy levels do not cross each 
other as a function of a for this system. This property, however, does not necessarily hold 
even for other two anyon systems and, in general, systems with three or more anyons in 
any potential have level crossings. 

Let us now combine the energy levels of the relative Hamiltonian with the energy levels 
of the CM system to obtain the full two particle energy levels - i.e., 



^particles = E CM + ^rel = (2j + p + 2 ± a/ 



(2.40) 
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Now compare these energy levels with the energy levels 

E n = (n + 1 V, n = integer (2-41) 

of a single particle in a two dimensional oscillator. For two particles, we would naively have 
expected the energy levels to be of the form 

E n ,m — (n+ 1)uj + (m + 1)uj n,m — integers. (2.42) 

This expectation is borne out only at a = and a = 7r, where it is clear that Eq.(2.40) is 
of the form of Eq.(2.42), so that two particle energy levels can always be written in terms 
of sums of single particle energy levels. However, for arbitrary a, the two anyon energy 
levels bear no simple additive or combinatoric relation to the levels of a one anyon system. 
This is the root cause of the difficulty in handling many anyon systems. 

A study of the ground state wavefunction also illuminates the same point. In the anyon 
gauge, the wavefunction is given by 

j>(R, 6, r, 9) OC e -™uj(R 2 +r*/4) r a/7r e i a e/ir = e -m^R 2 +r 2 /4) ^a/n ^ ^ 

where r in complex coordinates is x + iy = rcos9 + irsin9. This wavefunction can be written 
in the original two particle coordinate system as 

V>(ri, r 2 ) OC e -™>(rl+rl)/2( Ti _ ^a/* ^ ^ ^ 

which, incidentally, vanishes except for bosons (a — 0), thereby proving that anyons obey 
the Pauli principle. It is easy to see that for a = and a = ir, this wavefunction can be 
factorised into a product of two single particle wavefunctions — 

a = V(ri, r 2 ) = e — ?/2 e — 2 2 /2 = Xo(ri ) Xo(r2 ) (2 . 4 5) 

and 

^(ri,r 2 ) = rie -^?/2 e -^r 2 2 /2 _ -murl/2 e -murl/2 

a = n ' ' (2.46) 

= Xl(ri)xo(r2) - Xo(ri)xi(r2) 

where xo and xi are the ground and first excited states of the one particle system. For 

arbitrary a, however, the two particle wavefunction bears no simple relation to the one 

particle wavefunctions. Thus, many anyon states are not products of single anyon states. 

This is why even a system of free anyons needs to be tackled as an interacting theory. 
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There are other simple two anyon problems that can be solved exactly, such as the 
problem of two anyons in a Coulomb potential [15]. However, no three anyon problem 
has been completely solved so far. The main hurdle in solving systems with three or 
more anyons is that the relative phases between any two anyons depends also on the 
positions of all the other anyons in the system. This makes the problem quite intractable. 
Only partial solutions have been obtained so far using various different methods. Exact 
solutions for a part of the spectrum of three or more free anyons or in a harmonic oscillator 
potential (or equivalently an external magnetic field) have been found [16]. A semiclassical 
approach to the computation of the energy levels [17] has also been attempted. Perturbative 
approaches [18] in the limit of small a, where a is the statistics parameter, have also been 
used to demonstrate level crossings and the piecewise continuity of the ground state [19]. 
The level crossing phenomenon has been confirmed by recent numerical computations which 
have been used to obtain the first twenty odd energy levels [20] , which also showed that the 
analytic solutions found earlier formed a very small subset of the total number of solutions. 
Finally, there exist some exact results [21] regarding the symmetry of the spectrum as a 
function of the statistics parameter a. However, the computation of the full spectrum of 
energy levels of a three anyon system is beyond our capacity at the present time. Hence, we 
emphasize here that by and large, the study of three and more anyon quantum mechanics 
is still an open problem. 



2G 



Problems 

1. For a two-dimensional system of charged particles in a uniform magnetic field B 
(work in the gauge A = ^(—y,x) ), write down the Hamiltonian in terms of the 
complex coordinates z and z where z = (x + iy)/£ and z — (x — iy)/£ and £ is the 
magnetic length defined to be I = (hc/eB) 1 / 2 = {l/eB) 1/2 . 

a.) Show that [a, at] = 1 and [6, &t] = 1 where a = ^(2 + f ) and b = ^(2^ + f ). 

6.) Express the Hamiltonian in terms of the operators a,a\b and b* and find the eigen- 
values of H. What are the quantum numbers required to classify the states? How do 
you see the degeneracy of the states in terms of the quantum numbers? 

2. Solve the problem of Mott scattering of anyons (scattering from a 1/r potential). If 
you get stuck, look up Ref [15]. 
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3. Many Anyon Systems: Quantum Statistical Mechanics 



In this section, we shall study the quantum mechanical system of many anyons using 
the virial expansion of the equation of state, which is valid in the high temperature, low 
density limit. However, before starting on the many anyon problem, we shall show how 
the second virial coefficient of a system of free fermions or bosons is obtained and review 
in some detail, the cluster expansion method, which leads to the expression for the virial 
coefficients in terms of the cluster integrals [22]. 

For bosons and fermions, quantum statistical mechanics begins with the calculation of 
the grand partition function. The canonical partition function is given by 

Z(A,T) = J2 e ~ PEN (3- 1 ) 

states 

where En is the energy of the iV-particle system. In terms of single particle energy levels 
e p , En can be written as 

E N = ^2 e P n P ( 3 - 2 ) 

p 

where p enumerates the energy levels and n p is the occupation number of each level. 

The occupation number n p is just (0,1) for fermions and (0,1,2, ) for bosons, but it is 

constrained by the total number of particles given by 

N = J2n P - (3.3) 
p 

The canonical partition function cannot be easily evaluated because of this constraint. 
However, the grand canonical partition function defined by 

oo oo 

Z G (z,A,T) = J2 z N Z{A,T) = zN Yl e-0 E ' Cpnp (3.4) 

Spn p =JV 

where {n p } represents the collection of values of n p for the different states, invalidates 
the constraint, since it incorporates a sum over all N. Now, from Eq.(3.3), we see that 
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z N = U p z np and e ^ E p £ p n p = U p (e ^ £ p) n p yielding 

oo 

Z G (z,A,T) =J2Y, n(^ P ) np - ( 3 - 5 ) 

N=0{n p } P 

The whole purpose of considering the grand partition function instead of the partition 
function becomes clear when we realise that the two sums over N and {n p } reduce to 
independent sums over each n p -i.e., 

Z G (z, AT) = J2 \{{ze-^) n » = J] 5>e-*p)». (3.6) 

n p p p n 

Thus, the grand partition function is given by 

Z G (z, A, T) = + ze" 13 ^) for fermions, 

1 (3-7) 

and Zq(z,A,T) = „ — - for bosons, 

- LJ - 1 — ze~P e p) 
p v ' 

respectively. This method is inapplicable to anyons, because, as we have seen in Sec. (2), 
the energy levels of an iV-particle system cannot be computed in terms of the single particle 
energies. Hence, the exact calculation of the partition function for an anyon gas still remains 
an open problem. 

However, there exist other approximation schemes under which quantum gases are 
studied, one of which is called the cluster expansion (CE) method. The classical CE involves 
a systematic expansion of the interparticle potential, which is valid in the high temperature, 
low density regime. The quantum CE is defined by analogy with the classical CE. Using 
this expansion, corrections to the ideal gas law by systems of interacting quantum particles 
can be computed. This method is available to the anyon gas as well, because anyons can 
always be thought of as interacting bosons or fermions. 

Let us briefly review the cluster expansion method before applying it to the anyon gas. 
We consider a classical Boltzmann gas with interactions. The many particle Hamiltonian 
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is given by 

tf = £f^ + £^» r *-) ( 3 - 8 ) 

i i,j,... 

and the canonical partition function for N particles, obtained by integrating over all of 
classical phase space, is 

Z(A,T) = {2n) l Nm J d 2 n d 2 r N d 2 Pl <Pp N e-P H . (3.9) 

(Classically, the N particles are distinguishable and taken to be distinct. We divide by N\ 
to compensate for the overcounting.) The momentum integrations can be performed, since 
the integrand can be factorised into a product and the individual integrations are merely 
Gaussian -i.e., 

/d 2 p /3 P 2 ml .„ 
W? e ^ = w = i? (310 > 

where A is called the thermal wavelength. Thus, the canonical partition function, reduces 
to 

Z(A,T) = J d 2 r, d 2 r N e-^-V(^,-). (3.11) 

Now, the potential can be separated into sums of n-body potentials - i.e., 

V(ri, Tj, -) = E V(n, rj) + £ V{v h rj , r k ) + (3.12) 

We shall see that the lowest order correction to the ideal gas law involves only two body 
interactions, but let us keep the formulation general at this stage. The next step is to 
expand the integrand of Eq.(3.11) in powers of e^^ v — 1. We define 

, + (3.13) 

e -mr<**.) = (i + f ijk ) 

and so on. Each of these fijk... code for the n-particle interactions involved in the potential. 
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The partition function can be written in terms of the fijk.. as 

Z(A,T) = J d 2 n cPrNllil + h) ft (l + f ijk ) J] (1 • /;,/,/).... 

i<j i<j<k i<j<k<l 

= J d 2 n d 2 r N [i + (/ 12 + fa + ...) + (/12/23 + /12/14 + ...) + ...] 

+ (/l2/23/l3 + /l2/24/l4 + •••) + (/l23 + /l24 + •••) + ••• • 

(3.14) 

Now let us study, in more detail, each /-cluster involving interactions between / particles. 
A cluster integral b{ is defined as 

b\ = ^21-2^4 J ^ 2ri d 2 rN (contribution from /—clusters) (3.15) 

Explicitly, 



&! = - / d 2 Ti.l = 1, (3.1G) 



b 2 = TTvTT / ^ 2 ri^ 2 r 2 /l2 



2A 2 A 
1 

1 



d 2 ri d 2 r 3 f 13 (3.17) 



d Tid Tjfij, 



where the last equality is valid for any two arbitrary particles labelled i and j, and 

h = y d 2 vid 2 r 2 d 2 r 3 (f 12 fi3 + /13/23 + /12/23 + f 12/13/23 + /123) (3.18) 

with similar contributions from any other choice of three particles. Notice that the con- 
tributions obtained by changing the particle indices are all equal. Also note that the 
computation of 62 involves only two-body interactions even if the Hamiltonian includes 
three or more body interactions. 

It is now clear that the partition function may be written as a product of cluster 
integrals -i.e., 

Z(A,T) ~ U( nx2l ~ 2A ) mb ? 1 ( 3 - 19 ) 

{mi} I 

with YAmi = N. Here mi denotes the number of /-clusters and {mi} = (mi,TO2, ...) denotes 
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the collection of m/. However, we still need the combinatoric factor, which specifies the 
number of times a given cluster integral appears in the product. Let us compute this 
combinatoric factor. Firstly, since there are iV particles in the system, naively we expect 
iV! clusters. Hence, the R.H.S of Eq.(3.19) is multiplied by N\. However, particles in 
the same cluster are indistinguishable; so we have to divide by (l\) mi . Moreover, any two 
/-clusters are indistinguishable. Hence, for m/ /-clusters, we have to divide the R.H.S by 
to/!. Putting all this together, the partition function is given by 

W) = j^m E W^'-'^rj^—iT' 

{mi} I \ ■) i- 

1 Ah ( 3 - 2 °) 
' s X , mi\ A 2 

{mi\ I 

with Tlmi = N. (The simplest way, to convince yourself that this expression for the 
partition function in terms of the cluster integrals and all the combinatoric factors is correct, 
is to work out the partition function for a few particle systems explicitly.) As before, the 
constraint on the number of particles makes this partition function difficult to evaluate. 
Hence, just as was done for fermions and bosons, we move on to the grand partition function 
where the total number of particles N is also summed over. The expression for the grand 
partition function for this system is given by 



oo 

Z G (z,A,T) = J2^Z(A,T) 

N=0 



— 1 ( Abiz\ mi (3.21) 
m,P A 2 



n e ^(^) 

/ mi=0 

( Abyz 1 \ 



Notice that here again, the unrestricted sum over N and the sum over {mi} has been 
translated to an unrestricted sum over each of the m/. The equation of state obtained from 
this partition function is given by 
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Furthermore, the average number of particles < iV >= N is given by 

% lnZG = A 2 " 



N = z^\nZ G = ^Y,lhiz l . (3.23) 

l 



Combining Eq.(3.22) and Eq.(3.23), we get 



PA Zikz 1 biz + b 2 z 2 + ... 



NKT ^Ibiz 1 b lZ + 2b 2 z 2 + . 
Since b\ = 1, to first order in z, we find that 

PA 



(3.24) 



NKT 

Now, the virial coefficients of the system are defined by 

PA P 



l-b 2 z + 0(z z ). (3.25) 



-- Y aApX 2 ) 1 - 1 (3.26) 
NKT pKT 

where p is the density of particles. Thus, we have a power series expansion of the equation 
of state in terms of the inverse temperature A and the density p. This is clearly a useful 
concept at high temperatures and low densities, where only the first few terms in the series 
are likely to be relevant. Substituting for (pA 2 ) from Eq.(3.23), we see that 

oo 

pKT ^ lK l - x } (3.27) 

= ai + a 2 (z + 2b 2 z 2 + ) + ... 

to lowest order in z. Comparing this equation with Eq.(3.25), we see that a\ — 1 and 
a 2 = —b 2 . So to find the second virial coefficient, we only need to compute the two-cluster 
integral, which, in turn, depends only on two-body interactions. 

The quantum cluster expansion is defined by analogy with the classical cluster expan- 
sion. The quantum partition function is given by 

Z = ^2 rf2ri d 2 Y N i)* a (vi, ...r^e'^tpairi, ...r N ) 

a 1 (3-28) 
= mX 2 N J rf2 n d 2 r N W N (r h ....r N ), 

where {ip a } is a complete set of ortho normal wavefunctions for the system labelled by the 
quantum number a and the last equality defines Wn(ti, ...tn). By comparing this equation 
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with the definition of the classical partition function in Eq.(3.11), we see that 

W$(n, ...t n ) = e-fKMn*) = + /• •) (3.29) 

i<j 

(Since we shall only compute the second virial coefficient and we have seen that the second 
virial coefficient depends only on two body interactions, we have specialised to the case of 
two-body interactions alone.) Now, let us define new quantities Ui(r\, ...17) through 

W 1 (t 1 ) = U 1 (t 1 ), 
W 2 (n, r 2 ) = C/i(ri)C/i(r 2 ) + U 2 (r h r 2 ), 

(3.30) 

W 3 (ri, r 2 , r 3 ) = £/i(ri)[/i(r 2 )[/i(r 3 ) + £/i(ri)[/ 2 (r 2 , r 3 ) + C/i(r 2 )C/ 2 (n, r 3 ) 
+ L r i(r 3 )f/ 2 (ri, r 2 ) + C/ 3 (ri, r 2 , r 3 ), 

and so on. From Eq.(3.29), we see that the classical limits of the U\ functions can be 
identified as follows — 

Wf (ri) = 1 Uf( ri ) = 1, 
W 2 d (ri,r 2 ) = 1.1 + / 12 C/ 2 ( ri ,r 2 ) = / 12 , 
Wf ( ri , r 2 , r 3 ) = 1.1.1 + l.(/ 12 + / 13 + / 23 ) + (/ 12 / 13 + / 12 / 23 + / 13 / 23 ) + /12/23/13 
=>• [/^(ri, r 2 , r 3 ) = (/i 2 /i 3 + /i 2 / 23 + /i 3 / 23 ) + /i 2 /i 3 / 23 , 

(3.31) 

and so on - i.e., the t//(ri, ...r/) are the quantum analogs of the /-clusters in the classical case. 
Hence, for quantum statistical mechanics, the classical cluster integrals can be replaced by 
the quantum cluster integrals given by 

k = l[X 2l-2 A J rf2 n d2 rNUi(n, ...r,). (3.32) 

To calculate the second virial coefficient for any system, we need 6 2 , since we have 
already found that a 2 = — 6 2 , and to find 6 2 , we need to compute ^(ri,!^) which is a 
property of the two-body system. Let the Hamiltonian for the two-body system be 

H = -^l + ^l)+<\^-^\) (3.33) 
with eigenvalues E a and eigenfunctions ip a ( r i, r 2)- We transform to the CM (R) and 
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relative (r) coordinates to solve the problem. In terms of these coordinates, 

e iP-R 

^a( r l, r 2) — ► *Pa(R, r) = — y=-ip n (r) (3.34) 

and E a = P 2 /4m + e n . Here, the quantum number a has been split into (P, n) where P 
refers to the continuum quantum numbers of the CM system and n is the quantum number 
labelling the energy levels of the relative Hamiltonian. e n is found by solving the eigenvalue 
equation of the relative Hamiltonian given by 

[_I V 2 + u(r)]^ n (r) = e n (r)^n(r). (3.35) 
m 

The definition of W^(ri, r2) in Eq.(3.28) leads to its identification in this system as 
W 2 (r h r 2 ) = 2A 4 ]T ^ (n, r 2 ) e -^ a (n, r 2 ) 

2AA 4 /■ 9 e -0PV4m-/3e n 



/• -/3P 2 /4m-/3e„ 

y — i — i^wi 2 ( 3 - 36 ) 



(27T)2 

= 4A 2 $> n (r)| 2 e 
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which, in turn, identifies 

C/ 2 (ri,r 2 ) = 4A 2 ^ |^ n (r)| 2 e-^ - 1. (3.37) 

n 

Hence, the definition of the quantum cluster integral in Eq.(3.32) leads to 

b 2 = -a 2 = ^ / d 2 Rd 2 r[A\ 2 ^2\Mr)\ 2 e-^" - 1]. (3.38) 

The two individual terms in the integrand above give rise to area divergences when inte- 
grated over r and R. It is only their difference that is finite in the thermodynamic limit. 
(Even the difference is not always finite. The usual condition for finiteness of the differ- 
ence is that the interaction between particles should be short-ranged.) Often, it is more 
convenient to compute the second cluster integral as a difference between 6 2 for the system 
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under study and b 2 for some reference system, where b 2 is known exactly - e.g., b 2 could 
be computed when v(r) = 0. Then , for normalised relative wavefunctions -i.e., when 
f d 2 r\4> n (v)\ 2 = 1, we get 

b 2 - b° 2 = 2 J2( e ~ pen - e ~ Pel y ( 3 - 39 ) 



For the ideal Fermi and Bose gases, b 2 can be directly identified from the equation of 
state. For the Fermi gas, from Eq.(3.7), we see that the equation of state is given by 



PA 
~KT 



\nZ G = J rf2 plog(1 + 2e ^P 2 /2m ); {3AQ) 



since for free fermions e p = p 2 /2m. The R.H.S can now be expanded in a power series in 
z from which the b\ may be identified - i.e., 

PA _ A z 2 z 3 

A_. n; (3-41) 



i 

Hence, 62(fe rm io n s) is identified as 6°/ — —1/4- Similarly, for the Bose gas, 

PA A 



KT~ (2tt) 2 



J rf 2 plog(l-^ p2 / 2m ) 



A z 2 z 3 



(3.42) 



so that b 2 (bosons) is identified as b^ b = 1/4. 

We had earlier noted that the free anyon gas is already an interacting system, so a direct 
evaluation of the partition function is not possible. However, the second virial coefficient 
can be computed as long as we know the energy levels of the two anyon system. In Sec. (2), 
we explicitly computed the energy levels of two anyons in a harmonic oscillator potential. 
Let us use those results to find the second virial coefficient of an anyon gas in the same 
harmonic oscillator potential [11]. Ultimately, we shall take the limit where the oscillator 
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potential vanishes to obtain the virial coefficient of the free anyon gas. We shall choose our 
reference system to be that of free bosons (a — 0). From Eq.(2.39), we find that 



b 2 - b° 2b = 2j2[(j + i) e -W3+i+<*/*)u + je-Mj+l-aMu _ (j + 1)e -/?(2j+l)u, _ je -0(2j+l)u 
3=0 

oo 

= 4[cosh(a/7r - l)Pw - cosh j3u] je~ 2jl3uJ 

3=0 

cosh(a/7r — l)(3uo — cosh (3u 
sinh 2 (3u 

(3.43) 

To find the virial coefficient of a system of free anyons, we take the limit u — > 0, which 
yields 

*.-4 = k-) a --], (3-44) 

Z 71 71 



so that substituting b\ h = 1/4, we get 



6 2 = I[2(^-^ + l]. (3.45) 

4 7T 7T 



Notice that this equation is only valid for < a < 2tt, since a being a periodic variable 
specifying the statistics is periodic in 2n. Hence, the virial coefficient is non-analytic in a 
and has a cusp whenever a = 2nj for j an integer. By defining 5 such that a = 2nj + 5, 
we see that 

62 = j[2(-) 2 -- + l] (3.46) 

4 7T 7T 

is valid for any a. The virial coefficient can be plotted [10] as a function of a as shown in 
Fig.(13). 



37 




Fig. 13 



From the figure, it is clear that the values of the virial coefficient of an anyon gas interpolates 
between the values of that for fermions and bosons. 

Our result for the virial coefficient of an anyon gas is not obvious, though it is satisfying - 
just as spin and statistics of anyons are intermediate between spin and statistics of fermions 
and bosons, so are the virial coefficients. In fact, it is not even clear why the cluster 
expansion method works, because, when anyons are considered as interacting fermions or 
bosons, even the two-body interaction between particles is actually long range. However, 
the result that we have obtained here for the anyon gas is consistent, because the same 
answer has also been obtained using other regularisation schemes. Here, we used the 
harmonic oscillator potential as a regulator to discretise the energy levels and then took 
the limit where the oscillator potential vanishes to obtain the virial coefficient, whereas the 
original calculation of the virial coefficient involved a box normalisation [10]. 

The next logical step in this programme would be to compute the third virial coefficient. 
But, this would require the knowledge of three-body interactions. Since, as was mentioned 
in Sec. (2), the three anyon problem has not been completely solved in any potential so far, 
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the computation of the third virial coefficient too is an unsolved problem. However, as in 
the case of three anyon quantum mechanics, partial results have been obtained [23]. 



Problems 

I. In class, we studied the virial coefficient of a gas of anyons using a harmonic os- 
cillator potential as a regulator. Show that the same answer is obtained using box 
regularisation. (This was how it was originally done in Ref. [10].) 
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4. Many Anyon Systems : Mean Field Approach 



The basic idea of the mean field approach is to replace the effect of many particles by 
an 'average' or 'mean' field and to accomodate deviations from the mean field as residual 
short range interactions. In the context of anyons, the mean field approach involves replac- 
ing the flux-tubes carried by the charges by a uniform magnetic field with the same flux 
density [5] [6]. It is clear that this approximation is valid when the density of flux-tubes 
(equivalently particles) is high and fluctuations are small, -i.e., it is a high density, low 
temperature expansion. 

The many anyon Hamiltonian can be obtained by generalising the two anyon Hamilto- 
nian in Eq.(2.1) to N particles as 



N 



H = y y^~ g*) s (41) 



with 

Z X (Vj — V 



#3 



Thus the charge in each anyon sees the vector potential due to the flux-tubes in all the 
other anyons. We can also compute the magnetic field at the position of the i th charge. 
However, a naive computation leads to 

h = V x m = 0, (4.3) 

which is not surprising since a; can also be written as a gradient - i.e., 

where 9{j is the angle made by the vector (r, — rj) with an arbitrary axis. But using a 



40 



regularisation scheme [24] with 



a, = lim — 

e->0 27T f-f 



ZX {Ti- Yj 



Yi - Yj\ 2 + € 



2 ' 



(4.5) 



we can show that 



hi = V X a; = 0^5(Tj - Yj). 



(4.6) 



Not surprisingly, there is no magnetic field at the position of the i charge, due to any 
other particle, unless the two particles coincide. However, in the mean field approach, the 
flux-tubes are replaced by a constant magnetic field b with the same flux density. Let us 
assume that the density of anyons per unit area is given by p. Then the flux density of the 
anyons is given by 



when the integral is over a unit area and when p is sufficiently large. Hence, the appropriate 
uniform magnetic field to be used in the mean field approach is just 



so that we have a system of charges moving in a constant magnetic field as illustrated in 




(4.7) 



b = <f>p 




Fig.(14). 
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b ; at charges uniform b 
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Fig. 14 

Let us now study the quantum mechanical motion of a single particle in a constant 
magnetic field. The Hamiltonian is given by 

H = (P ~ q ^ (4.9) 
2m 

with the magnetic field B = V x A = constant. A can be chosen in many different ways 
(corresponding to different gauges), all of which lead to the same B. We shall work in the 
Landau gauge given by 

A x = -By, A y = 0. (4.10) 

(Another illustrative gauge (given as a problem in Chapter 2) is the symmetric gauge 
A x = —By/2, Ay = Bx/2, which is useful in the study of the Fractional Quantum Hall 
phenomenon.) Therefore, the appropriate Schrodinger equation that governs the motion 
of the particle is given by 

i iPx+ 2 q J y)2 + ^p^v) = E Px , P j Px , Py ^v)- (4.H) 

Here the p x and p y labels on the wavefunction and energy eigenvalues are the momentum 
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quantum numbers. Since there is no explicit x- dependence in the Hamiltonian, the motion 
in the x-direction is free and the wavefunction can be chosen to be of the form 

1>P.*y(*,v) = e ip * x X (y) (4.12) 

where x{v) satisfies Eq.(4.11) where, however, p x is now interpreted as the eigenvalue of the 
x-momentum operator. By defining qB/m = u (cu is the cyclotron frequency of a charged 
particle moving in the magnetic field B), and the magnetic length / = ^JTJqB = y/l/mw, 
Eq.(4.11) may be rewritten as 

\\vf + (f + P J) 2 ]x(y) = E Px , PyX (y)- (4-13) 

This is just the the Schrodinger equation of a shifted harmonic oscillator in the (y,p y ) 
co-ordinates. Hence, the energy eigenvalues are discrete and given by 

E Px , n = (n + l/2)u, (4.14) 

where we have replaced the continuum label p y in the subscript of E by the discrete numbers 
n. These discrete energy eigenvalues labelled by the integers n are called Landau levels. 
The corresponding eigenfunctions are given by 

X(y)= m 7== =e ^ H n (y/l + pJ) (4.15) 
vr ' V 2 n n\i 

where the H n are Hermite polynomials. Notice that the energy eigenvalues are independent 
of p x , which only effects a shift in the origin of the oscillator. Since for motion in a plane, 
p x is unrestricted, the Landau levels are infinitely degenerate. However, the degree of 
degeneracy becomes finite when the motion in the plane is restricted to a finite box with 
area A = L x L y . This degeneracy is easily computed. For motion in a one dimensional box 
of length L x , p x is quantised as 

p x = 2irn x /L x . (4.16) 
Furthermore, the allowed values of p x are restricted by the condition that the centre of the 
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oscillator has to lie between and L y . Hence, 

27m 



—^l 2 <L y (4.17) 



which, in turn, implies that the number of allowed values of n x , or equivalently, the degen- 
eracy of the Landau levels per unit area is given by 

n x 1 qB 



L x L y 2irP 2n 



(4.18) 



Let us now return to the many anyon problem which had earlier been reduced, in a 
mean field approach, to the problem of fermions or bosons moving in a uniform magnetic 
field. We shall choose our anyons to be fermions [6], rather than bosons, with attached 
flux-tubes, mainly because the problem of bosons in a magnetic field is itself unsolved. 
Hence, the analysis of fermions in a magnetic field (leading to anyon superconductivity) is 
easier. The other reason is that many properties of anyons appear to show a cusp in their 
behaviour in the bosonic limit. This was seen in Sec. (2), in the problem of two anyons in 
a harmonic oscillator potential, as well as in Sec. (3), where the second virial coefficient of 
the anyon gas was computed. Hence, the anyon gas may not have a smooth limit as the 
statistics parameter goes to zero if we start by perturbing from the bosonic end. However, 
the results that we shall obtain here by considering anyons to be fermion with attached 
flux-tubes have also been argued by starting with anyons as bosons [25] albeit with some 
approximations and hand-waving. Ultimately, whether we start with anyons as bosons or 
fermions is a matter of choice, and for each specific property of the anyon gas, one or the 
other approach may be easier. 

In the mean field approach, we need to solve the problem of fermions moving in a 
uniform magnetic field b = <pp = lapjq. (Remember that in Sec. (2), we had shown that 
the statistics factor a is given by q<f>/2.) The degeneracy of the Landau levels in this field 
is given by 

deg = £ = 5?. (4.19) 
Now, let us choose the statistics parameter a to be of the form a = n/n, where n is any 
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integer, so that the degeneracy is simply 

deg = ^. (4.20) 

n 

Since p is the density of particles and each level can contain p/n particles, clearly n Landau 
levels are completely filled. The next available single particle state is in the next Landau 
level which is an energetic distance uj = qb/m away. In the many-body or condensed matter 
parlance, this is called having a gap to single particle excitations. But if n is not an integer, 
then the last Landau level will only be partially filled and there will be no gap to single 
particle excitations. Hence, the parameter fractions a = ir/n appear to be special (e.g., 
like the magic numbers in the shell models of atomic and nuclear physics) and hence the 
states formed at these fractions should be particularly stable. 

To prove that the states formed at these special fractions are superconducting, we have 
to study the effect of adding a real magnetic field B to the fictitious magnetic field b and 
check whether a Meissner effect exists. The argument differs slightly depending on the 
relative signs of B and 6, and we shall consider both the cases separately. 

When the real magnetic field is aligned parallel to the fictitious magnetic field, they 
add and increase the degeneracy of the Landau level - i.e., 

deg = ?2±*>. (4.2!) 

But the number of particles per unit area p = nqb/2n remains unchanged. Hence, now the 
highest Landau level is only partially filled. Let us denote its filling fraction by (1 — x). 
From the conservation of density of particles, we have 

(n _ 1)S (^ + »(!+£) (1 _ l) = , = g (4 . 22) 

from which we see that 

(b + B)x = Bn. (4.23) 
Also from the energy eigenvalues in Eq.(4.14), we see that the total energy of p particles is 
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given by 

n-2 



„ g(& + -B)^,. L qib + B), 1 W1 , ,„ 



which can be simplified to give 



2tt m L 2 v 2 ; J y 1 

Substituting for x from Eq.(4.23), we get 

E = i^[^ + h -^-B\l--)} (4.26) 

For small external magnetic fields B, the energy relative to the ground state with no 
magnetic field is clearly positive and grows linearly with B. Thus, the anyon gas is a 
perfect diamagnet and tends to expel any external flux. 

When the external magnetic field B is in the opposite direction to statistical magnetic 
field b, the degeneracy of Landau levels decreases - i.e., 

deg = (4.27) 

So some of the p particles have to occupy the (n + l) th Landau level. Let us denote the 
filling fraction of the highest level by x. Then from conservation of particles, we have 

nq(b - B) q(b - B)x = nqb 

2vr 2tt 2tt 1 ' ' 

leading to 

(b - B)x = Bn. (4.29) 
The total energy of the system is given by 

E ='^-[tf + — -B\l + -)} (4.30) 

Notice that the linear term remains the same for B parallel and anti-parallel to b. Hence, 
once again for small B, the energy of this state relative to the state with B = is positive 
and grows linearly with B, emphasizing the need to repel any external magnetic field. 
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We have just demonstrated the Meissner effect by showing that the anyon gas finds it 
energetically favourable to exclude any external magnetic field. Thus, the anyon gas is a 
superconductor. The effect of adding an external magnetic field to the anyon gas can be 
depicted schematically as shown in Fig. (15). (In the figure, crosses denote particles and 
circles denote holes.) 
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Compare this with the schematic depiction (see Fig. (16)) of the creation of a quasiparticle- 
quasihole pair in the same system which costs an energy qb/m. 
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Fig. 16 

It is clear that the production of quasiparticle and quasiholes are closely related to 
the presence of a real magnetic field, because if the quasiparticle and quasihole were spa- 
tially separated, then the quasiparticle excitation is analogous to the situation depicted in 
Fig. (15c) with a real magnetic field —B and the quasihole excitation is analogous to the 
situation in Fig. (15b) with the field +B. (For p sufficiently large, the difference between 
p and p ± 1 is negligible.) With this connection, the Meissner effect can also be argued 
in the following way. Since for anyons, the magnetic field and hence degeneracy of Lan- 
dau levels is tied to the density, to accomodate any external magnetic field which changes 
the degeneracy, particles (or holes) have to be excited across the gap. This costs energy 
and hence, penetration by magnetic fields is unfavourable. Conversely, we also see that 
if particles do not fill Landau levels, there must be a real magnetic field to account for 
the mismatch between density of particles and degeneracy of Landau levels. Hence, this 
argument also implies that every quasiparticle excitation in the system is accompanied by 
a real magnetic field, so that in anyon superconductors, charged quasiparticle excitations 
and vortex excitations are indistinguishable, in contrast to usual superconductors where 



t 
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there are two types of excitations. 

To actually prove superconductivity in the anyon gas, we need to show that the col- 
lective excitation in the system is massless. This can be done by including fluctuations 
(residual interactions) about the mean field state [5] [6] in a random phase approximation. 
Relativistic field theory models have also been used to obtain the massless mode [26]. How- 
ever, even without going into the details of the calculation, there exists a simple heuris- 
tic argument to indicate the presence of the massless collective mode. Consider a very 
long wavelength density fluctuation (a collective excitation). Then p, though varying, is 
approximately constant over macroscopic lengths. Within each such macroscopic area, 
b = ap/n = p/n is constant. Hence, locally the system always has n filled Landau lev- 
els. Therefore, such an excitation neither requires any particle to be excited into a higher 
Landau level nor does it require any energy. Hence, such a wave is massless - i.e., we have 
proved the existence of a massless collective mode. 

Much further work has been accomplished in this field. Important questions like the 
persistence of superconducting currents at finite temperatures have been addressed [27], 
though the results are not yet conclusive. The connection between superconducting anyon 
states and FQHE states have been explored [28] and for interacting anyons, new super- 
conducting states have been found [29] at statistics parameter fractions a = n/u, where 
v is a fraction at which FQHE occurs. Also, the all-important question of experimental 
predictions and tests of anyon superconductors have been studied [30]. On the experimen- 
tal front, the one robust experimental prediction of all anyon models of superconductivity 
has been violation of the discrete symmetries parity P and time reversal T. But recent 
experimental results [31] (though somewhat controversial) appear to disfavour bulk P and 
T violation. Hence, the latest theoretical works [32] [33] have concentrated on studying 
the effects of layering on anyon superconductivity, in order to understand whether P and 
T violation survive layering. However, even if present day high T c superconductors are 
unlikely candidates of anyon superconductors, efforts to understand many anyon systems 
and, in particular, their superconducting behaviour remains undiminished. 
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Problems 

1. The standard mean field approach to anyon superconductivity starts from fermions 
in a magnetic field. But if we start with anyons as bosons plus attached flux-tubes, 
then the mean field problem to be solved is that of bosons (albeit hard-core because 
of the hard-core nature of anyons) in a magnetic field. Can you solve this problem? 
(This is as yet unsolved, but for starters, look up Ref . [34] . ) 
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5. Anyons in Field Theory 



Bosons and fermions are introduced in a second quantised field theory formalism 
through field operators that obey local commutation or anticommutation rules. However, 
anyon field operators , being representations of the braid group instead of the permutation 
group, cannot obey local commutation rules. This makes it hard to construct anyon field 
theories in the canonical way. Path integral quantisation is no easier, since the phase picked 
up by any anyon moving along a particular path depends on the position of all the other 
anyons in the system. Hence, unlike the case with bosons or fermions, paths cannot be 
weighted by a unique weight factor. But the study of anyons in a field theoretic formulation 
was rendered possible by the introduction [10] of an elegant idea called the Chern-Simons 
construction. Here, anyons were introduced as interacting fermions or bosons. It is this 
formulation (as opposed to more abstract formulations [35] ) that we shall study in this 
section. 

Let us consider any field theory in 2 + 1 dimensions described by a Lagrangian L and 
having a conserved current - i.e., d^j^ = 0. We can manufacture a gauge field a M and 
add to the Lagrangian 

AL = j^-^ a a^a a . (5.1) 

Is this an allowed extension of the Lagrangian? Firstly, if is to be a gauge field, AL has 
to be gauge invariant. Under — > — 

J AL d 3 x -> J ALd 3 x- J j^K) d 3 x + ^ J e^ Q (^A)(^a Q ) d 3 x. (5.2) 

The third term obviously vanishes (upto surface terms) on integrating by parts due to 
the anti-symmetry of e^ va . The second term also certainly vanishes on integrating by 
parts at the equation of motion level, since d^j^ = 0. However, for any explicit current, 
we can, in fact, construct a gauge invariant coupling of the current with a local gauge 
field. For example, when is a fermionic current, the gauge invariant coupling is just 
3n a ^ — ^l^^- But when is a scalar current, j^a^ is replaced by D^tpD^ip. Hence, the 
AL in Eq.(5.1) is certainly an allowed gauge invariant extension of any Lagrangian. 
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The term (AL) C s = (/i/2)e^ a a^a Q is called the Chern-Simons (CS) term. The CS 
term and its non-abelian generalisations are interesting field theories in their own right 
and have been studied for years [36] in other contexts. More recently, they have shot 
into prominence as prototypes of 'topological' or 'metric independent' field theories [37]. 
However, for our purposes here, we shall only recollect the following salient features of 
Chern-Simons theories. An abelian CS field theory is described by 

L = —F^F^ + ^ a A»d»A a , (5.3) 

where F 0t = cfoA — diAo = E{ and (l/2)eoyi™ = B. It is clear that the electric field E{ 
is a two component vector and the magnetic field B is a pseudoscalar, and that these are 
the only non-zero components of F^, since we are in 2 + 1 dimensions. The parameter 
ji in Eq.(5.3) has the dimensions of a mass and is the gauge invariant mass term for the 
gauge field. This can be seen by explicitly computing the propagator. Furthermore, the CS 
term is odd under the discrete symmetries parity P (x — > —x, y — > y) and time reversal T 
(t — > —t). Thus, the Lagrangian in Eq.(5.3) describes a U(l) gauge theory with a massive 
photon. A gauge invariant mass for the photon has been introduced at the expense of the 
violation of P and T. 

With this introduction to CS* or topological field theories, let us get back to the study 
of Eq.(5.1). Notice that the Lagrangian in Eq.(5.1) does not include the usual kinetic piece, 
the F tli) F lxv term for the gauge field. Also, from the equation of motion, we get 

Jn = fie^ad"^, (5.4) 

which, when coupled with a gauge condition allows for a solution for the gauge field in terms 
of the current j^. ( We shall see this explicitly later in this section.) Hence, the motion 
of the gauge field is completely determined by the current and has no independent 
dynamics. In this respect, CS theories of relevance to anyons (or metric-independent field 
theories) differ from most of the CS'-models studied earlier. However, the gauge field does 
affect the statistics of the current carrying particles. Integrating the zeroth component of 
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Eq.(5.4) over all space, we get 



J eoy(dV) 



(5.5) 



q = fuj). 



Thus, every charge q of the current is accompanied by a flux and is an anyon. This 
mechanism of attaching fluxes to charges is called the CS construction. The current j M 
can be a fermion number current (j^ = ^7^) m which case fermions turn into anyons, 
or a bosonic current (j^ = (j^d^cj) — (d^ycj)) so that bosons turn into anyons, or even a 
topological current so that topological objects like solitons and vortices turn into anyons. 



Let us compute the statistics of these particles in the field theory context. In the path 
integral formalism, when two such objects (flux carrying charges) are exchanged, the phase 
is given by e lScyi where S ex is the action involved in exchanging them. 
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Now, it is clear that S ex = (l/2)S TO t, where SVot is the action for taking particle B all 
around particle A. 
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Fig. 18 

S TO t, in turn, is twice the Aharanov-Bohm action Sab required to take a single charge 
around a single flux, because S ro t involves taking both a charge and a flux around a charge 
and a flux. Therefore, 

S ex = Srot/2 = Sab- (5.6) 

The Aharanov-Bohm phase is easily computed since we know the 'electromagnetic' La- 
grangian (Eq.(5.1) ) that governs the motion of a charged particle in the field of a flux-tube 
and is given by 

S cx = Sab = / dt 

r u r (5 - 7) 

= / jia 1 dt+ - eoijdidodj dt 

(in the ao = gauge). Since ji — Q v i, the first term in Eq.(5.7) gives 

J jia 1 dt = q J a • dl = q<j>. (5.8) 

Also, from the equation of motion in Eq.(5.4), we see that the motion of the gauge field is 
related to the current as 

{ie 0ij d a J = -ji. (5.9) 
Substituting Eq.(5.9) in Eq.(5.7), we see that 

S ex = q(P/2 = a. (5.10) 
Thus, the statistics of a charge with a flux induced on it by a CS construction is a = q<j)/2, 
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consistent with what we had earlier derived for the charge-flux composite in the Hamiltonian 
formulation. 

The relation between the Lagrangian CS field theory formulation and the Hamiltonian 
formulation for many anyons introduced in Sec. (4), can be explicitly demonstrated [6] by 
carrying out the canonical transformation between the Lagrangian and Hamiltonian for- 
mulations for a gas of anyons. We start with the Lagrangian formulation of a gas of anyons 
(represented by point particle bosonic charges with attached flux-tubes) given by 

L = Y^^ + ^o(r«) + gr Q -a(r Q )]- | J e^d v a a d 2 v (5.11) 

where a is the particle index and the coupling between the particle with charge q and 
the gauge field a M is as in standard electromagnetism. However, the usual kinetic term of 
electromagnetism F^ V F^ V is now replaced by the CS term, in accordance with Eq.(5.1). 
This Lagrangian may be rewritten as 

L = ^2^-i- 2 a + J aoO'o- neoijdidj) d 2 r + g^r Q - a+ | J eoijaidj d 2 r (5.12) 



where 



J0 = g^5 2 (r-r a ) (5.13) 



represents the point particle charge density of the anyons. Notice that apart from the first 
term, all the terms are linear in ao or time derivatives. Also, since there is no kinetic term 
for ao, the equation of motion with respect to ao yields the constraint 

dL 

= => j = fieoijdiaj = fib, (5.14) 

so that the field strength is completely determined by jo- But for an abelian gauge 
field theory, the entire gauge invariant content of the gauge field a% is contained in the field 
strength fy. Hence, if we eliminate the extra degrees of freedom in a,{ by a gauge condition, 
we can actually solve for aj in terms of jo, so that the motion of the gauge field is entirely 
determined by the fields forming jo and has no independent dynamics. The CS gauge field 
has been introduced merely to attach flux-tubes to charges and turn them into anyons. 
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To go to the Hamiltonian formulation, we require the canonical momenta given by 

dL 

mr a + q& (5.15) 



rut n • 

and 

dL pi 

dtii 2 

The Hamiltonian can now be written as 



Pi — wt- = —'— I e 0ij a J d 2 r. (5.16) 



H = ^ a + PiCti - L 
a 

= ^(mr a + qa) ■ r a - | J eoijina.j d 2 r - ^ -j r l ~ Q XI *" ' a + f J e 0y d 



2 r 



m ■ 2 (Pa - ga) 



2 



a a 

(5.17) 

which, in terms of the velocity is just the Hamiltonian of free particles. So the classical 
equations of motion are identical to the equations of motion for free particles. What 
has been altered, however, is the relation between canonical velocity and momenta (see 
Eq.(5.16) ). Hence, the quantum commutation relations are no longer the same. This is 
the significance of introducing the CS term without the usual kinetic piece and is consistent 
with our original introduction of in Sec. (4), as just a way of enforcing anyon boundary 
conditions in a different gauge. 

Let us now solve for aj in terms of jo from Eq.(5.14) and the additional gauge condition 

<9;a* = 0. (5.18) 
The solution to these two equations, for the jo given in Eq.(5.13), is given by 

a * (r) = 2^ l^'^y^f^') 



(5.19) 



which can be easily checked. The gauge condition is obviously satisfied due to the anti- 
symmetry of eoij. Also, by regularising the denominator of aj as explained in Sec. (4), we 
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can explicitly check that 

b = V x a = / j'o(r / )27r5(r - r') dV = (5.20) 

2-Kfl J fl 

Notice that the many body Hamiltonian in Eq.(5.17), along with the solution for a« in 
Eq.(5.20), is precisely the many body Hamiltonian (Eq.(4.1)) and gauge field a,{ (Eq.(4.2)) 
that was used in Sec. (4), where it was obtained by generalising the two anyon Hamiltonian. 
Thus, we have established the equivalence of the Lagrangian CS formulation and the many 
body Hamiltonian formulation used in the earlier sections [6]. 

Let us now study a study a specific example of a field theory - the abelian Higgs 
model with a CS term [12] - whose solitons (classical solutions of the equations of motion) 
are 'anyonic'. We shall first construct the topologically non-trivial vortex solutions of the 
abelian Higgs model and then show that they are charged and have fractional spin in the 
presence of a CS term. The Lagrangian for this model is given by 

L = -\F^ + \{d^- iqAJ<l>*{&> + iqA»)<P - c 4 (0*0 - ^) 2 + ^ a F^A a . (5.21) 

This model has a U(l) gauge symmetry. However, when C2 > 0, the potential energy V(4>) 
is minimised when 

V{4>) = c 4 (0*0 - ^) 2 = =► |0| = ^ = v, (5.22) 

where v is the vacuum expectation value. Hence, the U(l) symmetry is spontaneously 
broken by the vacuum. 

The usual vacuum has 

0(r, 6) = v, A(r, 6) = and A (r,9)=0, (5.23) 

which is a solution of the equations of motion and perturbation theory is built up by 
expanding the fields in modes around this vacuum. But besides the vacuum solution, this 
theory also possesses topologically non-trivial finite energy solutions [38]. To have finite 
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energy solutions, all we need to ascertain is that 



\<j>(r -> oo,0)| = v, A(r -> oo,0) = and A (r -> oo, 0) = 0, (5.24) 

so that the energy integral does not diverge. However, the condition for the scalar field 
is satisfied even if the scalar field has a non-trivial phase at infinity. So a solution of the 
equations of motion which satisfies 

<j>(r ^oo,9)= ve in6 

Id ■ n n (5.25) 
=^d 8 <l>ar^oo,e))=~ve tne = i- 

r ad r 

has finite potential energy. Finiteness of the kinetic energy of the scalar field also implies 
that 

J (di - iqA i )(j)*(d i - iqA { )<\) d 2 r = finite, (5.26) 
which, in turn, gives the conditions on the asymptotic behaviour of A and Aq as 

Tl 

A e (r -> oo,6) = — , A r = and Aq = 0. (5.27) 
qr 

Solutions that satisfy Eq.(5.25) and Eq.(5.26) are topologically non-trivial vortices. They 
are called vortex solutions because they carry flux - i.e., using these solutions, we see that 
the flux is given by 

J B rf 2 x = J A • dl = J {^)rd6 = (5.28) 

where n is called the vorticity. 

These solutions are topologically stable because solutions with different values of n are 
not deformable into each other without changing the scalar field configuration throughout 
the infinite boundary of space. This argument is easily understood pictorially. In Fig. (19), 
we have the vacuum solution with = v everywhere, including at spatial infinity, 



58 



t tttt 



Fig. 19 

whereas in Fig. (20), we have the one-vortex solution with 0(oo,#) = ve 





Fig. 20 



Both the vacuum and the one- vortex configurations satisfy V(4>) = 0. But to deform one 
solution to another, we need to go through configurations where 0(oo, 6) is neither v nor 
ve iG and which have V(4>) ^ 0. Since this has to occur throughout the infinite boundary 
of two dimensional space, this will cost an infinite energy. In general, from Eq.(5.21), we 
see that the vacuum expectation value v takes values on a circle (since is complex and 
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2 = v 2 ). Hence, the configuration space of the vacuum is denoted by S^. The boundary 
of real (2-dimensional) space is also a circle and is denoted by S^. Hence, the boundary 
condition expressed in Eq.(5.25) denotes a mapping from to and the vorticity n 
represents the number of times encircles when r encircles S}. once. Each of these 
solutions falls in a different topological class and is stable, since it requires an infinite 
energy to change the configuration of <fi throughout the infinite boundary. 

Now, let us introduce a, CS term for the gauge field. We had earlier seen that in 
the presence of the CS term, the flux and charge get related. Hence, now, these vortex 
solutions also possess a charge given by 

077-71 

Q = rf> = v(— -). (5.29) 

Since these charged vortices carry both charge and flux, by our earlier arguments in Sec.(l), 
they are anyons with spin j = Qcp/Air = im 2 n/q 2 and statistics phase a = Q(p/2 = 
2ir 2 n 2 /q 2 . Notice that here an explicit kinetic term for the gauge field has been included, 
so that the gauge field is really a dynamical degree of freedom. Hence, this model is not 
really relevant to anyons. However, even when the usual kinetic piece is switched off, the 
model continues to exhibit anyonic solutions [39]. 

This was just one example of a field theory whose solitonic excitations are anyons. 
Another example is the 0(3) c-model with the Hopf term, whose skyrmionic excitations 
are anyons [28]. In all such models, the basic ingredient is the CS term (or equivalently 
the Hopf term) which relates the charge and the flux. 

The CS term can also be extended to non-abelian theories with the appropriate La- 
grangian being given by [36] 

L = UiF^ - ^ a tr(F^A a - 2 -gA»A v A a ). (5.30) 

This Lagrangian can be shown to be gauge invariant provided 4iifi/g 2 = integer. As in 
the abelian theory, \x is the gauge invariant mass for the gauge boson and the CS term 
is odd under parity and time reversal. The pure non-abelian CS Lagrangian without the 
usual kinetic piece has been recently related to problems in topology and knot theory, 
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integrable models in statistical mechanics, conformal field theories in 1+1 dimensions and 
2+1 dimensional quantum gravity. (For references to the original papers, see Ref. [28].) 
The remarkable feature of pure CS theories that makes it amenable to exact solutions 
is its general covariance, — the Lagrangian is not only Lorentz invariant, but it is also 
generally covariant without any metric insertions. Hence, correlation functions of the pure 
CS theory depend only on the topology of the manifold and not on details such as the 
metric on the manifold. Besides theoretical interest in non-abelian theories, there have 
also been speculations [40] that non-abelions, — generalisations of anyons that form non- 
abelian representations of the braid group, — may play a role in the even denominator 
FQHE. However, since non-abelian CS theories are not directly related to anyons, we shall 
not pursue this topic any further here. 
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Problems 

1. Similar to vortices and charged vortices in 2 + 1 dimensions, in 3 + 1 dimensions we 
have monopoles and dyons. The Georgi-Glashow model in 3 + 1 dimensions, which 
consists of an 5*0(3) gauge field A® interacting with an isovector Higgs field is one 
of the simplest examples. The Lagrangian for this model is given by 

L = —trF^F^ + ±tr DrfD^ - - a 2 ) 2 , 

where = d^A a v - d„A£ - gt ahc A\A c v and {D^) a = d^(f) a - ge ahc A b ^ c . Can you 
construct monopole solutions for this field theory by analogy with the vortex solution 
in Eqs. (5.25) and (5.26)? 

b. Why are monopoles topologically stable? 

c. Can you guess the generalisation of the Chern-Simons term that induces charge on 
the monopoles (converting them to dyons)? 

Some useful references are Refs. [41] and [42] . 
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6. Anyons in the Fractional Quantum Hall Effect 



Any study of anyons would be incomplete without an account of its most outstanding 
success — its application to the FQHE. This is the only physical system where there exists 
incontrovertible evidence for the existence of anyons, because quasi-particle excitations over 
the FQHE ground state have been explicitly shown to obey fractional statistics. However, 
since familiarity with the Quantum Hall system is not in the repertoire of the average 
graduate students, we shall introduce the background material in some detail. 

Let us first remind ourselves of the Hall Effect. Here, electrons in a plane show trans- 
verse conductivity when a magnetic field is applied perpendicular to the plane. The Hall 
geometry is depicted in Fig. (21). 
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Electrons are allowed to move from a source S to a drain D causing a current J, which is 
measured, as are the longitudinal and transverse voltage drops Vl and Vjj- The existence 
of a non-zero Vjj can be explained just by classical electrodynamics using the equation 



F = eE + ev x B 



(6.1) 



for the electrons with charge q — e. For the geometry in Fig. (21), we have the equations 

eE T eVvBn . cv x Bq 



whose solution is 



• Jx ev y c>Q 

V X = + — 

m m 



v x = voe tuJt , v v = -— 1 + ivoe 



, v y 



m 



iujt 



5, 



o 



(6.2) 



(6.3) 



where vq is the initial velocity and u = eBo/m is the frequency of the cyclotron motion. 
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The constant term in v y represents the drift velocity and the Hall current is given by 

peE x 

JH = pev d j Hx = 0,j Hy = — , (6.4) 

-DO 

where p is the density of charge carriers. By changing the Hall geometry, we can also find 

3Hx = ^ET 1 , 3Hy = 0, (6.5) 

-DO 

when the electric field is along the y-direction. Hence, if we define a conductivity matrix 

Oij by 

/ rr rr \ 

(6.6) 

then classically, a xy = —a yx = —pe/Bo and a xx = a yy = 0. The resistivity matrix 
Pij = a^j 1 and Eq.(6.6) leads to 

Pxy = -Pyx = — and p xx = p yy = 0. (6.7) 
ep 

We have not taken collisions between electrons into account in deriving the above classical 
result. However, a more realistic picture requires collisions and by incorporating them, we 
get the semiclassical equations of motion given by 



/ ®xx 


®xy 


Oij = 1 




V Oyx 


a yy 


—pe/Bo an d 


®xx 



. eE x q < v y > B < v x > 
<v x > = 1 



and < Vn, > 



mm r 
e < v x > bo < v y > 



(6.8) 



y 

m 



where < v > is the average velocity of the electrons and r is the average time between 
collisions. In equilibrium, < v >= and we get 



re/m ujr 2 e/m 

( w ¥ + i) Vy = ~(u^ + iy 



v x = , -,x E x, v y = - j — E x (6.9) 



where, as before, u = eBo/m. Note that unlike Eq.(6.5) which does not give the correct 
Bq — > limit, (Hall current goes to infinity instead of vanishing), here u — * as Bq — > and 
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the Hall current vanishes. Thus, in the semiclassical limit, the elements of the conductivity- 
matrix are 

pe 2 r/m pe 2 r 2 uj/m 



(wV + l) ^ ° yX ~ ° Xy - (cu 2 T 2 + l) 



and the elements of the resistivity matrix are 



(6.10) 



m muj Bq 

Pxx = Pyy = o and Pyx = p X y = 7) — • 
pe z r pe £ ep 



(6.11) 



Notice that the off-diagonal elements of the resisitivity are unchanged from their classical 
values. 

Now, when the Hall resistance was experimentally measured [43] in 1980, at low tem- 
peratures (0-2 degK) and in strong magnetic fields (1-20 Tesla ), a surprising result was 
found, an showed plateaux (as shown schematically in Fig.22) instead of varying linearly 
with 1/Bq as expected classically or semiclassically. 
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Furthermore, the longitudinal resistance, instead of being constant as expected by Eq.(6.10), 
vanished at the plateaux and peaked in between as shown (schematically) in Fig. 23. 
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Even more surprisingly, though the behaviour of the Hall resistance at the edges of the steps 
was non-universal, the Rjj at the plateaux was unexpectedly constant and reproducible to 
the accuracy of one part in 10 5 . This effect was called the Integer Quantum Hall Effect 
(IQHE), because the midpoints of the plateaux occured at integer values of a ratio v = pj Pb 
(where pb = gBq/2t[). Since pb is the degeneracy of each Landau level, this ratio which 
expressed the fraction of filled Landau levels was called the filling factor. The filling factor 
could either be thought of as a measure of -Bo at fixed p or as a measure of p at fixed Bq. 

The remarkable accuracy of the experimentally measured resistances is explained by 
the following observation. Usually the resistances Rh = Vjj/I and Rl = Vl/I are related 
to the resistivities by geometric factors of length or width. However, for the Hall geometry, 
when L » W and the voltage drops are measured sufficiently far from the actual edges 
so that the applied current density is uniform, we see that 

RH = — = ^ W =P X y. (6.12) 

Hence, geometric factors which can never be measured to accuracies of one in 10 5 have 
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cancelled out leaving the transverse resistance equal to the transverse resistivity. In fact, 
the quantisation of the experimentally measured Hall resistance is so accurate, that the 
Quantum Hall system is now used for the most precise determination of the fine structure 
constant. 

The midpoint values of the plateaux can be easily understood by studying the quantum 
mechanical problem of non-interacting electrons in a transverse magnetic field. In Sec. (4), 
we studied the problem of an electron in a transverse magnetic field and found that the 
energy levels are given by 

eBo 



E n = {n +1/2)- 



m 



with the degeneracy 



Deg = p B 



2tt ' 

where we have substituted u = eBo/m in Eqs.(4.14) and (4.18). Hence, when the filling 
factor u = p/ pb is an integer, it is clear that precisely an integer number of Landau levels 
are filled. As mentioned in Sec. (4), this implies a gap to single particle excitations, given 
by uj = eBo/m. Also, unlike the case in Sec. (4), here Bq is independent of p and there 
is no argument for a massless collective excitation. On the contrary, explicit calculations 
show that the collective excitation is massive. So the system is particularly stable when the 
density (or equivalently the magnetic field) is such that v is an integer. From Eq.(6.11), 
we see that for these values of the density, 

B Bp 1 2tt 1 

Rh = — = — = , 6.13 

ep epB integer e z integer 

in accordance with the experimental values of Rh at the midpoint of the plateaux. Thus, 
the densities of the electrons at the midpoint of the plateaux have been identified with fully 
filled Landau levels and the correct Hall conductivity is predicted for these densities. 

However, to understand why the conductivity remains fixed even when Bq (or p) is 
changed from the midpoint value is harder. At a hand-waving level, we can argue that 
because of the stability of the system when v is an integer, even when the field Bq is 
changed slightly, the system prefers to keep the average density fixed such that v is an 
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integer and accomodate the deviation in Bo as a local fluctuation. These local fluctuations 
do not contribute to the conductivity because they get 'pinned' or 'localised' by impurities 
in the sample. Hence, the conductivity stays fixed at the value that it had for v = integer. 
At a slightly more rigorous level, the idea is that weak disorder in the system (due to 
impurities and imperfections in the sample) leads to the formation of some localised states 
, whereas other states are extended. Current can only be carried by extended states. Hence, 
if the density of states as a function of the energy had the pattern shown in Fig. (24), 
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Fig. 24 

then as the Fermi level spans each localised region, the current will remain constant (plateau 
in the Hall conductivity), while when it spans the extended states, the current (and hence, 
conductivity) will increase. To study in more detail the kinds of disorder potentials that 
can lead to the density of states diagram in Fig. (24) would lead us too far afield. For 
further details, we refer the reader to Ref.[44]. 

In 1984, in stronger magnetic fields and cleaner samples, the FQHE was seen [2], where 
the midpoints of the plateaux were found to occur at fractional values of the filling factor 
such as 1/3, 1/5, etc. This was unexpected because the single particle quantum mechanics 
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analysis that was used for the IQHE predicted that at v — fraction, the system would 
be highly degenerate and not at all stable. So, the very existence of the effect showed 
that inter-electron interactions which were ignored in the earlier study must be important. 
Hence, the appropriate Hamiltonian for the FQHE is expected to be 

i i i^j 

where V(rj — r-,-) denotes the inter-electron Coulomb repulsion and V(rj) is a neutralising 
background potential. Unlike the IQHE problem, this Hamiltonian cannot be solved exactly 
even in the absence of disorder. 

There have been several approaches to the theoretical understanding of the FQHE, 
mainly because no theory has provided a complete picture yet. Here, we shall concentrate 
on two approaches, both of which appear to use anyon ideas in a fundamental way — the 
trial wavefunction approach and the Landau-Ginzburg-Chern-Simons field theory approach. 

The trial wavefunction approach was initiated by Laughlin [3] who proposed that the 
ground state of the FQHE system is described by the wavefunction given by 

N 

Mz!, Z2, ....z N ) = Y[( Zi - Zj re-^ 2 / 412 (6.15) 

i<j 

for the fraction v = 1/m. Here z% = x% + iyi is the position of the i th particle and m has 
to be an odd integer to satisfy the criterion that the particles are electrons. He arrived at 
this wavefunction from certain general principles such as a.) the wavefunction should be 
antisymmetric under exchange, b.) the wavefunction should comprise of single particle states 
in the lowest Landau level, c.) the wavefunction should be an eigenstate of total angular 
momentum, and d.) inspired guesswork. This ansatz for the ground state wavefunction has 
had enormous success, mainly because no other ansatz has been found with lower energy. 

To understand the significance of this wavefunction, let us write 

H 2 = e"^ (6.16) 
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so that 




(6.17) 



By setting the fictitious temperature (5 = 1/m, we see that <f> can be interpreted as the 
potential energy of a two dimensional gas of classical particles of charge W repelling each 
other through a logarithmic interaction and being attracted to the origin by a uniform 
(opposite) charge density pu = l/2nl 2 . This potential energy is minimised (i.e., \tp\ 2 is 
maximised) by a uniform distribution of the charge 'm' particles and is 'electrically' neutral 
everywhere, when the average density of the charge 'm' particles is precisely equal to pjj, 
which in turn implies that the average density of the electrons p is equal to l/2nml 2 . 
Hence, at these densities, the many electron wavefunction peaks when the coordinates Z{ 
are uniformly distributed and is expected to be energetically favourable, since the Coulomb 
repulsion in Eq.(6.14) is minimised. For other values of p, the classical gas has excess charge 
'm' particles either near the origin or near the sample boundary and hence fails to be 
uniformly distributed. So, the appropriate ip peaks at a configuration of the z^'s that does 
not describe a uniform distribution of electrons and suffers from a high repulsive Coulomb 
energy. Hence, the Laughlin picture is that as the density is reduced from pb = eB/2n 
where one Landau level is filled, whenever p = p m = ps/m =>■ v — 1/m, the energy is 
minimised and the state is stable. 

The Laughlin state at v = 1/m can be shown to have massive fractionally charged 
quasiparticle and quasihole single particle excitations, which obey anyon statistics. To see 
this, let us insert adiabatically one unit of flux through an imaginary thin solenoid piercing 
the many body state. The intermediate Hamiltonian as the flux is being evolved through 
the solenoid is given by 

»=£ (p, ~ e tT A ' )2 + E v ' (r '» + \ £ v ' (r - - r -» (fU8) 

i i i=£j 

where is the gauge potential of the flux through the solenoid. As the flux changes, the 
wavefunction also changes so as to be an eigenstate of the Hamiltonian. At the end of the 
process, however, the Hamiltonian has returned to the original Hamiltonian, since one unit 
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of flux through an infinitely thin solenoid can be gauged away. (Remember that the term 
eA^ in the Hamiltonian can be traded for a phase e te< ^ in the wavefunction in the 'anyon' 
gauge. But when the flux <fi = <J)q = 2n/e, the phase is just 1 and hence irrelevant.) But 
the state has not returned to the original state. As flux is adiabatically added, every single 
particle state 

z m e -\z\ 2 /4 z m+l e -\zf/4 (og) 

-i.e., its angular momentum increases by one. The state with the highest angular momen- 
tum moves over to the next Landau level and a new state appears at m — 0. If a single 
unit of flux is removed adiabatically, then 

z m e -\z\ 2 /4 ^ z m-l e -\z\ 2 /4_ (a20) 

Hence, a state from the next Landau level moves down and the m = state disappears. 
Thus, the effect of adding (or removing) one quantum of flux and then gauge transforming 
is to increase (or decrease) the angular momentum of the single particle states by one unit. 
But if we assume that the original state described by the Laughlin wavefunction is non- 
degenerate, (since no other states with the same energy have been found), then the new 
state, after evolution of the flux, has to describe an excited state of the original Hamilto- 
nian, with a higher energy eigenvalue. This automatically proves that the quasiparticle or 
quasihole excitations have a gap. 

The electric charge of these excitations can also be easily computed. Let us assume 
that the flux has been evolved through a very thin solenoid. Faraway from the solenoid, 
we expect that this state will be indistinguishable from the ground state, except that every 
level of the single particle states has moved over to the next level. Hence, if the flux point 
is surrounded by a large circle, then the charge that has entered or left the circle is just 
the average charge per state, provided the total charge and the total number of states are 
uniformly distributed. In position space, single particle states are labelled by the position 
vector r and are uniformly distributed with a degeneracy of eB/2n. By the plasma analogy, 
charge is also uniformly distributed in real space, with the density ep m = e 2 B/2irm for 
v — 1/m. Thus, the charge per state is just ±e/m which is identified as the charge of the 
excitation. 



71 



The statistics of these quasiparticles and quasiholes can be explicitly found [45] by a 
Berry phase calculation. However, a much simpler way is to notice that these excitations 
have a flux 27r/e and a charge e/m and hence are anyons with fractional spin and statistics 
given by 

q<ft elm x 2-n/e 1 qcp elm x 27i7e n . . 

3 — i— — : = t; — an d ot — — — = — . (6.21) 

4tt 4tt 2m 2 2 m v 1 

Eq.(6.19) implies that an explicit ansatz for the wavefunction of a quasihole excitation can 
be written as 

V^o, *i, = 11^ - 11^ - ^) m e" E<N2/4/2 (6.22) 

where is the position of the infinitely thin solenoid. This ansatz is expected to be good 
except very near the solenoid. The wavefunction for the quasielectron involves derivative 
operators and is not as straighforward to understand, at least within the Laughlin scheme. 

Laughlin's wavefunctions only explained the plateaux at the fractions v — 1/m where 
m was an odd integer. But, experimentally, plateaux were seen at many other rational 
fractions v = p/q where q was an odd integer. To explain the other fractions, the hierarchy 
scheme was evolved [4]. The idea was that quasiparticle or quasihole excitations over the 
Laughlin state themselves behaved like particles in a magnetic field and could form new 
correlated many body states which could represent the FQHE state at other fractions. 

The hierarchy scheme is easily understood in the anyon language. Firstly, notice that 
the wavefunction for a quasihole excitation remains analytic (see Eq.(6.22)). Also, quasi- 
holes obey anyonic statistics. Hence, the simplest possible wavefunction for a collection 
of many quasiholes of generic charge qe and statistics a' (using Laughlin's arguments to 
arrive at his wavefunction, but now remembering that the particles are anyons) is given by 

M 

^(z 01 ,zo2,....^oM) = n(^-^) 2fc+Q 'e^ (|g,/e)El|2wi2/4Z2 . (6.23) 

i<j 

This wavefunction looks exactly like Laughlin's wavefunction except that m is replaced 
by 2k + a' to account for the changed statistics and in the exponent I 2 is replaced by 
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el 2 /\q\ to account for the charge of the quasiholes which is qe. The same plasma analogy 
now suffices to find the density of quasiholes p q h for which this wavefunction describes a 
uniformly distributed electrically neutral plasma and is hence likely to be an energetically 
favoured wavefunction. The plasma has 'charge' (2k + a') particles repelling each other 
and being attracted to the origin by a uniform 'charge' density p = \q\/2irel 2 . So by our 
earlier argument, we shall have a uniform density of quasiholes when 

= W^) x 2S' (6 - 24) 

Now, we know that the charge of the quasiholes is qe = — e/m and their statistics parameter 
a = n/m =^> a' = 1/m. (Remember that we are using complex notation in Eq.(6.23). 
Hence, under exchange of particles, (z re \) a = (re t0 ) a — > ( re l ( fil + 7r )) a — (z m \) a e ma ). 
Using this, the filling fraction of the quasiholes is given by 

v qh = ^ = l —r- (6-25) 

PB m{2 k + —) 
m 

To find the equivalent density of electrons, notice that at a given density of electrons, the 
total charge remains fixed, whether it is counted as quasihole charges or electron charges. 
Hence, the total charge carried by the quasiholes is given by 



qp q h = -%qh = — — — (6-26) 



m 



which, in turn, is equal to ep eq ^, where p^ is the equivalent density of electrons. So the 
total density of electrons and hence, the filling fraction is given by 

v = ^ l ~~r = —^r- ^ 

m 2 (2k-\ ) m H — - 

m 2k 

For m = 3 and k — 1, v — 2/7 and for m = 3 and k — 2, v — 4/13. Thus, we have obtained 
other fractions with odd denominators. This process can now be iterated. We can consider 
excitations over the v = 2/7 or the v = 4/13 state and form new correlated states with 
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those excitations. A little algebra shows that by repeating this procedure, we get filling 
fractions which can be written as 

v= ^ . (6.28) 

m H 



2fci + - 



2k 2 + 



1 



2k s 

Quasielectron excitations have a statistics factor 2k — a instead of 2k + a and have the 
same sign of charge as the electrons. Repeating the same analysis as above for quasielectron 
excitations gives v as in Eq.(6.28) except that all the plus signs are replaced by minus signs. 
In general, when we allow for both quasiparticle and quasihole excitations over each state, 
the possible filling fractions are given by 

v = ^ . (6.29) 

m + W) 

2k x + - 2 



2k 2 + ...+ aS 



2k s 

where ai are either +1 or -1 depending on whether quasiparticle or quasihole excitations 
are involved. All rational fractions with odd denominators are obtained once in this way. 
Also for the FQHE system, the hierarchy, as the above scheme to generate the fractions 
is called, works in the sense that for any fraction that has been observed, all the other 
fractions that lie before it in the hierarchy have also been observed. 

The problem with the hierarchy scheme is that some fractions (e.g., v = 5/13, at the 
third level of hierarchy) have not been observed. Also, for other fractions like v = 6/13 
which is seen at the fifth level of hierarchy starting from the v — 1/3 state, the number of 
quasiparticles of various types is so much more than the number of original electrons, that 
it is hard to understand how the explanation of the original electrons forming a uniformly 
distributed v — 1/3 state survives. A way out of this predicament was suggested by 
Jain [13], who could directly obtain the wavefunctions for all the odd denominator fractions 
that have been seen. The starting point of his approach was to note that the FQHE 
was phenomenologically very similar to the IQHE, and hence, the theories of both the 
phenomena should also be related. 
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Let us start with an IQHE state, say, the p-filled Landau level state represented by 

„ = ^=p =► B = ^, (6.30) 
PB pe 

and then attach fluxtubes of strength = 2/c0o to each electron. We know that this 
will convert fermions to anyons, but provided k is an integer, the statistical parameter 
a = e0/2 = 2kn is irrelevant (we have chosen the statistical charge to be the same as the 
real electric charge, which is always possible) and the electrons remain fermions. These 
fermions with attached fluxtubes are called 'composite fermions'. Now, within the mean 
field approach (valid for a high density of fermions), the flux of the fluxtubes can be spread 
out and we have ordinary fermions moving in an effective magnetic field 

2-ko 4nk 

B eS = ±B + 2k<p = ±— F - + , 6.31 

pe e 

so that the filling factor v for ordinary electrons is given by 

i = P_B^SM3l = ± l +2k= 2 j!E±l. (6 .32) 
v p p p p 

Thus, the basic idea is that the FQHE for ordinary electrons occurs because of the IQHE for 
composite electrons — i.e., the same kind of correlations between electrons are responsible 
both for the IQHE and the FQHE. When p — 1, this procedure yields the Laughlin fractions 
v = l/(2k± 1). 

This idea can also be used to write down trial wavefunctions for arbitrary odd denom- 
inator fractions. For instance, to obtain the Laughlin wavefunctions, we start with the 
wavefunction for one filled Landau level given by 

...z N ) = Jl(zi - Zj)e-^\ 2 /* 12 (6.33) 

i<j 

and then adiabatically introduce fluxtubes with <fi = 2k<po at the site of each electron. We 
know that adiabatic evolution of a flux unit at any point zq involves the factor Hi(zi — zq) 
from Eq.(6.22). Hence, insertion of 2k fluxtubes at the positions of all the electrons leads 
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to the wavefunction given by 



ikk+i(zi, ■■■■zn) = JJ( 

i<j 



Zi-Zj) ipi(zi,...z N ) 



(6.34) 




i<j 



which is precisely the Laughlin wavefunction for the filling fractions v = l/(2k + 1). The 
same procedure can also be used to write down wavefunctions for the other fractions v = 
p/ (2kp + 1) as 



where if) p (zi, ...zn, z\, ■■■zn) is the wavefunction for p filled Landau levels. But the form 
of the wavefunctions at v — p/(2kp — 1) is not as obvious. Explicit quasiparticle and 
quasihole wavefunctions can also be written down very easily in this formalism. Just as the 
FQHE wavefunctions are obtained from the IQHE wavefunctions ip p , the quasiparticle or 
quasihole wavefunctions for the FQHE are obtained from the quasihole wavefunctions or 
quasielectron wavefunctions ip~ of the IQHE, by multiplying them by the factor Ui < j(zi — 
Zj) 2k . The quasihole wavefunction for the v — 1/ (2k + 1) state in this formalism, coincides 
with the Laughlin quasihole for the same state, but the quasielectron wavefunction involves 
higher Landau levels and differs from Laughlin's ansatz and other trial wavefunctions in the 
literature. However, it appears much more natural and less arbitrary and works quite well 
numerically. The quasiparticle and quasihole charges and statistics can also be computed as 
was done for Laughlin's wavefunctions and the hierarchy wavefunctions and they coincide 
with the earlier results for any fraction v — p/q. 

The wavefunction approach is a microscopic approach and has been fairly successful 
in explaining the phenomenon of the FQHE. However, it is not completely satisfactory, 
because it fails to illuminate all the symmetries and does not provide a complete un- 
derstanding of the problem. To give an analogy, it is as if soon after superconductivity 
was discovered, the iV-body projected BCS wavefunction in the coordinate representa- 
tion was directly written down. Although it was correct, a complete understanding of the 



Z i -Z j ) 2 %(z h ...ZN,Z 1 ,...Z N ) 



(6.35) 



i<j 
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phenomenon of superconductivity would not have been possible without discovering the 
Landau- Ginzburg theory and the phenomenon of Cooper pairing which actually led to the 
BCS theory Hence, for the FQHE system too, there have been several attempts to find 
analogues of the Cooper pair that condenses and a Landau- Ginzburg theory. It is only in 
the last couple of years that these attempts have begun to bear fruit and an effective field 
theory called the Landau- Ginzburg-Chern- Simons (LGCS) theory has been evolved [14]. 
This approach uses anyon ideas and the Chern-Simons construction introduced in Sec. (5) 
to write down an effective action for the FQHE problem, whose saddle point solutions cor- 
respond to FQHE states. The key to this approach is the realisation that within the mean 
field approximation, real magnetic field is indistinguishable from the fictitious or statisti- 
cal magnetic field introduced by a Chern-Simons term. Hence, real magnetic flux can be 
moved on to fermions and identified as statistical flux, so that the fermions get converted 
to anyons. The interesting result that was discovered in this approach was that precisely 
at those values of the magnetic field where the FQHE occured, all of the real magnetic 
flux could be transferred onto the fermions so that they could be transmuted to bosons. 
Pictorially, this procedure can be represented as shown in Fig. (25). 




fermions in a 
Magnetic Field 



fermions + flux tubes 



Fig. 25 



free bosons 



An equivalent picture (which is the one that we shall implement formally) is to consider 
the fermions as bosons with attached flux tubes. Then the FQHE occurs precisely when 
the external magnetic field cancels (in a mean field sense) the effect of the flux tubes. Once 
again, this can be depicted as shown in Fig. (26). 
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fermions in a 

Magnetic Field fermions + flux tubes free bosons 

in a Magnetic Field 

Fig. 26 



In either picture, FQHE occurs when the fermions turn into free bosons on the average. 
Thus, the FQHE (or formation of an incompressible fluid) is equivalent to the formation of 
a Bose condensate. The stability of the FQHE states is 'explained' by the well-known fact 
that bosons can lower their energy by Bose condensing. 

Let us now derive the above explanation more formally, starting from the microscopic 
Hamiltonian for electrons in an external electromagnetic potential given by 

N N 

H = ^J2&- eA ( r *)) 2 + E eA ^) + E v ^ - r ^ ( 6 - 36 ) 

i=l i=l i<j 

where V{ji — rj) is the repulsive interelectron potential. (In writing this Hamiltonian, we 
have dropped the background neutralising potential, which, of course, is always present.) 
The solution ■0(ri,r2, ...r/v) to the Schrodinger equation 

H^(r h r 2 , ...rjv) = E^(n, r 2 , ...r N ) (6.37) 

has to be antisymmetric under the exchange of any two coordinates, since the particles are 
fermions. These fermions can also be written as bosons interacting with an appropriate 
Chern-Simons gauge field. Hence, an equivalent formulation of the problem is given by 

N N 

H ' = ^ E(p* - eA ^ - ea *) 2 + E + E v ^ - r ^)' ( 6 - 38 ) 

i=l i=l i<j 

with 

H'<f>(n, r 2 , ...r N ) = E<j>(n, r 2 , ...r N ). (6.39) 
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(f> is now symmetric under the exchange of coordinates, and sl{ in Eq.(6.38) is given by 



where #y is the angle made by (rj — rj) with an arbitrary axis. 

To prove that the two systems are really equivalent, all we need to show is that the factor 
a /en in Eq.(6.40) is precisely the factor needed to convert fermions to bosons. Consider a 
unitary transformation on the wavefunction and the Hamiltonian in Eq.(6.39) — i.e., 

0(ri, r 2 , ...r N ) — ► 0(n, r 2 , ...rjy) = C/0(ri, r 2 , ... rjV ) = e _ ^ Ei <^0(ri, r 2 , ...rjy) (6.41) 



and 

if' — ► UH'U- 1 = H. (6.42) 
The transformed Schrodinger equation is given by 



UH'U^U^n, r 2 , ... rjv ) = EU<j>(n, r 2 , ... rjV ) 
=>• H<p(r 1 , r 2 , ...rjv) = £0(ri, r 2 , ...r^). 



(6.43) 



Furthermore, since 9{j is the angle made by (rj — rj) with an arbitrary axis, 

Bij = Op + 7T. (6.44) 

Hence, if (f> is symmetric under exchange, picks up the phase e ta and is fermionic whenever 
a = (2k + l)n. For these values of a, the systems described by Eqs.(6.36) and (6.37) on the 
one hand and Eqs.(6.38) and (6.39) on the other, are equivalent — i.e., 0(ri, r 2 , ...r/v) = 
■?/>(ri,r 2 , ".rjv). This is not really surprising, since the unitary transformation in Eq.(6.41) 
is just the generalisation to many particles of the gauge transformation used in Sec. (2), to 
go from the fermion gauge to the anyon gauge, a = (2k + l)ir is the value of the statistics 
parameter for fermions to become bosons. 
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We can also explicitly check that the Hamiltonian in Eq.(6.38) can be obtained from a 
Chern-Simons Lagrangian. From Eq.(6.40), we can show that 



b( Ti ) =Vxa=^ 8(ri - Tj) = (6.45) 



where p is the density of particles. Comparing with Eq.(5.5), we see that this is precisely 
the form of the equations of motion derived from a Chern-Simons Lagrangian. Hence, the 
appropriate action that describes the effective field theory of the FQHE problem is given 
by 

S — S a + Sfj, 
r e 2 

where S a = d 3 x{ — e^a^ap} 

and Sa= I d 3 x{^[id t - e(A + a )U + n$<j> + ^[^(-iV - eA - ea) 2 ^]} 
J 2m 

d 2 nd 2 r 2 {0t( ri )0t( r2 )y( ri _ r 2 )0( ri )0(r 2 )}. 

(6.46) 

This field theory formalism now enables us to look for minima of the action (which corre- 
spond to the usual mean field theories) and incorporate quantum corrections by expanding 
around the minima. 

Let us first consider the case where there exists a non-zero magnetic field, but no electric 
field so that 

A = 0, and e lj diAj = -B = constant. (6.47) 
Here, the action in Eq.(6.46) is minimised by choosing 

(j) — y/p, A = - < a >, and a = 0, (6.48) 



where the density of particles p is a constant, {(jr = p is enforced by choosing the chemical 
potential p suitably). Since the statistical magnetic field b is related to p (Eq.(6.45)) and 
Eq.(6.48) relates b to the external magnetic field Bq, this minimisation is only possible 
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when 



or when 



B Q = b = ^ (6.49) 



v= — = 9 =-= (Q 50) 

V p B eBo/2n a (2k + 1) 1 ' ' 

where ps = eB/2n = degeneracy of each Landau level and we have substituted a = 
(2k + Thus, the action is minimised at the densities for which the filling fraction 

v = 1/odd integer -i.e., the Laughlin fractions. To prove that the vacua at these fractions 
are incompressible, we also need to show that all excitations over these vacua are massive. 
For the quasiparticle excitations, we have already seen in Sec. (5) that the CSLG theory 
or the abelian Higgs model with a CS term has charged vortex solutions that have finite 
(non-zero) energies and fractional spin. Collective excitations, which are fluctuations of 
(A + a) (which could be massless, in principle, as for anyon superconductivity) are also 
massive because of the spontaneous breakdown of the U(l) symmetry caused by the vacuum 
expectation value of 0. Hence, it appears reasonable to identify the classical minima of the 
effective field theory with the Laughlin states. 

Let us now calculate the Hall conductance by applying an external electric field E{ = 
— d{Ao along with the external magnetic field Bo = —e l;1 diAj. The observed current can be 
obtained from the action by using 

. _ ds _ds^ _ds^ 

since A{ and enter symmetrically. The equations of motion with respect to en are 
given by 

OS dS a dSs 

— = 0^-^ + -± = 0. (6.52) 
Substituting for dS^/ddi in Eq.(6.51), we see that 



e 2 

Ji = 



Now, from the minimum energy ansatz in Eq.(6.48) for a static magnetic field and the 
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ansatz E{ = diA$ for the external electric field, we have 



doaj = —doAj = and <9jao = —djAo = Ej. 



(6.54) 



Thus, the current can be written as 



Ji 



2ir(2k + l) 




(6.55) 



so that we obtain the Hall conductance as 



xy — 



2n(2k + 1) 



(6.56) 



which agrees with the semiclassical answer for the Hall conductivity a xy = ep/B given in 
Eq.(6.5), since p/B = ev/2-n and v = I /(2k + 1). Hence, the LGCS theory 'explains' the 
stability of the FQHE states at the Laughlin fractions and gives the right Hall conductivity 
at these densities. 

The charm of the LGCS theory lies in the fact that besides reproducing the phe- 
nomenology of the FQHE, it also provides a formalism for addressing questions like the 
existence of an order parameter and off-diagonal long range order in the system. The full 
implications of the LGCS theory are yet to be understood, although the theory has been 
taken considerably further [46]. This is an active field of research and many more questions 
remain, both to be formulated and answered. 



1. One of the points which was glossed over in the lecture was the question of how 
states get 'pinned' or 'localised' by impurities. Just to get a feel for this point, solve 
the quantum mechanical problem of a single fermion in a magnetic field and in the 
presence of an impurity scattering potential Vj = \5(x — xo)5(y — yo). (If you get 
stuck, look up Ref.[47].) 

2. A subtlety. 

Fermions with attached flux tubes of strength 2/c0o remain fermions because the 



Problems 
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statistical parameter a = (§)0 = (§)2/c(-^) = 2kn is defined only modulo 2ir. So 
why are 'composite fermions' different from ordinary fermions and lead to different 
physics (FQHE vs. IQHE)? 



83 



7. Conclusion 



We conclude this review by pointing out that these lectures are merely meant to serve 
as an eye-opener to the ever-expanding field of anyon physics. Besides, the several 'known' 
unanswered questions in the field, there probably remain many more unexplored and un- 
expected connections between CS theories and other topological and non-topological field 
theories. Applications of anyon ideas to other phenomena in condensed matter physics also 
remain a distinct possibility. Our hope is to inspire many more readers to join the 'anyon 
bandwagon'. 
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